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Abstract 

We present rigorous results for several variants of the Hubbard model in the strong- 
coupling regime. We establish a mathematically controlled perturbation expansion 
which shows how previously proposed effective interactions are, in fact, leading-order 
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terms of well defined (volume-independent) unitarily equivalent interactions. In addi- 
tion, in the very asymmetric (Falicov-Kimball) regime, we are able to apply recently 
developed phase-diagram technology (quantum Pirogov-Sinai theory) to conclude that 
the zero-temperature phase diagrams obtained for the leading classical part remain 
valid, except for thin excluded regions and small deformations, for the full-fledged 
quantum interaction at zero or small temperature. Moreover, the phase diagram is sta- 
ble upon the addition of arbitrary, but sufficiently small, further quantum terms that 
do not break the ground-states symmetries. This generalizes and unifies a number of 
previous results on the subject; in particular published results on the zero-temperature 
phase diagram of the Falikov-Kimball model (with and without magnetic flux) are ex- 
tended to small temperatures and/or small ionic hopping. We give explicit expressions 
for the first few orders, in the hopping amplitude, of these equivalent interactions, and 
we describe the resulting phase diagram. Our approach, however, yields algorithms to 
compute arbitrary high orders. 
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1 Introduction 

1.1 Scope of the paper 

In this paper we present rigorous results for several variants of the Hubbard model in the 
strong-coupling regime. For a general class of Hubbard-type models, we construct effective 
Hamiltonians with explicit exchange interaction terms by using unitary conjugations ob- 
tained from convergent perturbation expansions. Furthermore, we determine ground states 
and low-temperature phases in the asymmetric (Falicov-Kimball) regime. Our results are 
applications of tools recently developed in [§, || for the rigorous study of quantum statis- 
tical mechanical lattice systems. These tools can be applied to systems with finite-volume 
Hamiltonians of the form H\ = ifoA + V\, where: (1) the dominant part, Hq\, is "classical" 
—in the sense of being diagonal in a tensor-product basis (see Section below) — and has 
a known energy spectrum with spectral gaps uniformly positive in the volume A C 2^, and 
(2) the perturbation, Va, must be small in a sense to be made precise below, but can be 
rather general. In particular, it may involve interactions of infinite range as long as some 
exponential-decay condition [ (|2.20|) below] is satisfied. For the (one-band) Hubbard-type 



models to be studied in this paper, ifoA involves only on-site terms, and the perturbation is 
just nearest-neighbor hopping. These models rank among the simplest applications of our 
conjugation method. To illustrate a slightly more complicated application, we also analyze 
the 3-band Hubbard model. 
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The methodology used in this paper has two parts: (I) a rigorously controlled perturbation 
expansion, the "conjugation method", adapted to the statistical mechanical setting, and (II) 
a theory of stability of phases and phase diagrams known as quantum Pirogov-Sinai theory. 
Part (I) yields a perturbative expansion with mathematically controlled error terms, which 
agrees, to leading orders, with well known expansions |5], ^7], |39], ||, [7|, [L4|, |31[j. Part 



(II) allows us to extend — and rederive in a more systematic fashion — rigorous results on 



low-temperature phase diagrams for models such as the Falicov-Kimball model |26|, IIR |16 



30, 25 



In the statistical mechanics of classical lattice systems, there exist efficient methods and 
well tested intuitions to determine and describe ground states and equilibrium states. The 
mathematical study of quantum lattice systems is less advanced. The loss of the probabilistic 
framework — noncommuting observables, complex-valued expectations — necessarily makes 
an analysis more abstract, and often there is a divorce between the simple and concrete 
intuitions practitioners resort to and a mathematical approach suitable for rigorous proofs. 

The study of quantum perturbations of classical systems offers an excellent opportunity 
to establish a bridge between the classical and quantum formalisms and to extend techniques 
and intuitions from the classical realm to the quantum realm in a mathematically controlled 
manner. This is, in essence, what our methods accomplish. 

Of course, the scope of such methods is limited. They are inadequate to rigorously study 
strong quantum-mechanical correlations, such as those appearing in Fermi liquids or super- 
conductors, or to analyze the spontaneous breaking of continuous symmetries. To be more 
specific, our perturbation expansions and quantum Pirogov Sinai theory converge only if the 
quantum perturbation V\ is small, in a sense explained later in this paper, as compared to 
the classical Hamiltonian H ^. But the perturbation expansion enables us to decompose the 
Hamiltonian if a into an effective classical Hamiltonian and a quantum perturbation in a way 
that takes into account how the original perturbation V\ lifts degeneracies in the energy spec- 
trum of the original classical Hamiltonian Hq\. In particular, it enables us to show how — as 
a consequence of the Pauli principle — effective exchange interactions are generated, starting 
from a Hamiltonian with spin-independent interactions, such as the Hubbard Hamiltonian. 
The key idea underlying our perturbation technique is to perturbatively construct a unitary 
conjugation of the original Hamiltonian which block-diagonalizes it up to some finite order in 
the quantum perturbation. It is analogous to Nekhoroshev's method in classical mechanics. 
The goal of the method is to unitarily conjugate the original Hamiltonian to an effective 
Hamiltonian of a form that enables us to apply quantum Pirogov-Sinai theory. As a result, 
we are able to study genuine quantum effects. 

The main purpose of this paper is to illustrate the perturbation technique (I) and quan- 
tum Pirogov Sinai theory (II) by analyzing concrete models, such as Hubbard-type models, 
of interest in condensed matter physics. Further applications can be found in |S], |T2|] . 

We emphasize that the two methods, (I) and (II) above, used in this paper are indepen- 
dent of each other. The perturbation method can be applied whenever the leading classical 
Hamiltonian has an appropriate, essentially volume-independent discrete low-energy spec- 
trum and its ground states simultaneously minimize all local contributions to the classical 
Hamiltonian (technically, the classical Hamiltonian is assumed to be defined by a so-called 
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m-potentials [21] — see Section |2|), and the quantum perturbation is "small" (see Section g). 
It may or may not happen that the transformed Hamiltonian, the "effective" Hamiltonian, 
is in a form suitable for the application of quantum Pirogov Sinai theory. If it is, then we 
are able to rigorously investigate the low-temperature phase diagram in certain regions of 
the space of thermodynamic parameters. But even if it is not in a form enabling us to apply 
method (II), it may nevertheless yield new insight into the structure of ground states or 
low-temperature equilibrium states (by enabling us to e.g. appeal to "educated guesses"). 
As mentioned, the effective Hamiltonian may, for example, include exchange interactions 
between spins - absent in the original Hamiltonian - that enable one to determine, at least 
heuristically, the type of magnetic ordering of the ground states. We are confident that our 
perturbation method will be useful in providing suitable effective Hamiltonians as starting 
points for alternative approaches to the study of phase diagrams, such as renormalization 
group methods. 

A number of delicate concepts are involved when designing mathematically controlled 
approaches. Some of them are discussed below (Sections |2j] and |4~1~1 ). We would like to 
start here with a general remark related to the infinite- volume limit: As we are after statis- 
tical mechanical properties (phases, phase transitions), we have to pass, eventually, to the 
thermodynamic limit. This implies that we need to cope with Hamiltonians for arbitrarily 
large volumes A, e.g. arbitrarily large squares. To achieve this, we shall adopt the following 
two "principles": 

(PI) We shall attempt to establish estimates or bounds that hold uniformly for sufficiently 
large volumes. 

(P2) Instead of working at the level of Hamiltonians, we shall work with so-called interac- 
tions, which are families of local operators {$>x :IcF finite} (so-called |A|-body 
terms), in terms of which the Hamiltonians are given by if a = Y^xnk^^x- Interac- 
tions are defined without reference to the total volume of the system. Hence operations 
performed on interactions automatically apply to all finite- volume Hamiltonians. 

Next, we proceed to introducing some of the models studied in this paper and to sum- 
marizing our main results. 



1.2 Models and results 
1.2.1 One-band Hubbard model 

We shall study the one-band Hubbard model and some of its variants. These are models 
of two types of fermions — which, for concreteness, will be called spin-up and spin-down 
electrons, but which can also be considered to be different species of particles, like ions and 
spin-polarized electrons — whose quantum dynamics, on a finite square lattice A, is governed 
by the Hamiltonian 



Hk = U n x+ n x - n x+ - fi_ J] n x _ 
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"EE It 1 ?" 1 City* + flT 1 . (1.1) 



Here c|. CT and are electron creation- and annihilation operators satisfying the usual an- 
ticommutation relations [see (|2.9|) - (|2.10|) 1. The number operator of an electron of spin 



a = +, — at the site x is n xa = c xa c xa . The chemical potentials of the electrons with 
up-spin and down-spin are denoted by fi + and //_, respectively. The symbol (xy) (resp. [xy]) 
denotes an unordered (resp. ordered) pair of nearest neighbor sites of the lattice. 

The Hubbard model is the simplest model believed to embody the essential physics 
governing strongly correlated electron liquids. It exhibits the interplay between Coulomb 
repulsion and the kinetic energy term of electrons. The model is widely studied in the context 
of phenomena such as magnetism and high-T c superconductivity, which are believed to result 
from such an interplay. It is a lattice model of interacting fermions in the tight-binding 
approximation. The model retains only that part of the (screened) Coulomb repulsion which 
manifests itself when two electrons of opposite spin occupy the same lattice site, the strength 
of this interaction being given by the coupling constant U. The kinetic energy term describes 
the hopping of electrons between nearest-neighbor sites. We allow a direction-dependent 
hopping: the hopping amplitude of electrons of spin a from a site x to a site y is denoted by 
the symbol t^ x \ For the Hamiltonian to be self-adjoint, t^ must be the complex-conjugate 
of t^ . Complex hopping amplitudes are encountered in the presence of an external magnetic 
field. 

If = t 1 ^ the Hamiltonian ( |1.1| ) does not change under a rotation of the spin quan- 
tization axis. Hence the model has an SU(2) symmetry. This is the case for the standard 
one-band Hubbard model for which t+ = = t, for all nearest neighbor sites x, y. The 



highly asymmetric regime t+ <C will be called the Falicov-Kimball regime. The 

limiting case = constitutes the Falicov-Kimball model (with magnetic flux |Jl8|). In 
this setting, up-spin electrons are interpreted as static ions, while down-spin electrons are 
viewed as scalar quantum-mechanical particles. If the magnetic flux is zero, £_ is real and 
independent of x, y. 

We study the model ( [Lip in the strong coupling limit U ^> . This corresponds 
to the situation in which the Coulomb repulsion is much larger than the bandwidth of 
(uncorrelated) electrons. The large on-site repulsion forbids double occupancy of a site at 
zero temperature. In this limit, the terms in the first line on the RHS of Ql.l|) act as a 
leading classical part, H$\, while the hopping term can be treated as a perturbation, V\. 
As we shall see, this gives rise to a relatively simple application of our perturbation scheme, 
due to the on-site character of the leading interaction. 

To discuss our results, we first present the zero-temperature phase diagram of the dom- 
inant part, Hqa, of the Hamiltonian, in the plane of the chemical potentials. This is a 
purely classical interaction, and ground states are tensor products of single-site spin states 
minimizing each on-site term. It is given in Figure |I| The symbols {+}, {— }, {±} and 
{0} are used to denote ground states of iJ A in which each lattice site is occupied by an 
up-spin electron, a down-spin electron, two electrons of opposite spins and no electron, re- 
spectively. All the boundary lines in this phase diagram are lines of infinite degeneracy; 
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e.g. for < jj + = /j- < U the ground state is infinitely degenerate, because all singly 
occupied configurations are equally likely. The point /i + = \i_ = (orig in) and the point 
/x+ = /i_ = U also correspond to infinitely many ground states. At the origin, each site is 
either empty or singly occupied, whereas, at the point (U,U), each site is either singly- or 
doubly occupied. Next, we describe our results for the Hubbard model. 



(I) Controlled perturbation expansion 

We restrict our attention to those values of the chemical potentials for which the spectrum 
of the Hamiltonian H \ can be decomposed into two spectral bands separated by an energy 
gap which is large compared to the hopping amplitudes. We discuss mainly the half-filled 
regime (Section ||), for which the lower band corresponds to the subspace of states with singly 
occupied sites. In this regime, to have a large-enough spectral gap, we must consider a region 
in the positive quadrant in the plane of chemical potentials in which /i < fi- < U — /io, 
where /io 3> tjf^ . This is the shaded region in the plane of chemical potentials shown in 
Figure 0. 

For values of the chemical potentials in this shaded region, we determine, for each finite 
order n, an interaction j that is 

(Rl) equivalent to the one corresponding to ( p..l|) , i.e. it gives rise to Hamiltonians related 
to (|1 . lp by a unitary transformation, and 



in the sense that the 



(R2) block- diagonal to order t n /U n , where t = sup^ 

matrix elements of the operators $^ between the lower and higher bands (of Hq\) are 
of the prescribed order. 

Transformations of this sort have been devised before, see for instance || |20], |27|, [39], [| 
[7], [TJ], 31], though not in a manner suited for rigorous analysis, as we comment upon below, 



in Section \I7l\ Our expressions, obtained with the technique developed in 0, yield, when 
restricted to the lower band, the well known leading order terms: Heisenberg interaction 
+ terms of order t A /U 3 favoring valence-bond states (and breaking the degeneracy among 
those) + terms of order t 6 /U 5 + ... In particular, they coincide exactly with those of [ j3lfl , 
in the half-filled regime. (Our approach can be viewed as a mathematical justification of 
the series proposed there.) Nevertheless, each transformed interaction {^^j has terms, 

^x 7^ 0' with arbitrarily large |X| but correspondingly small order in t/U, which have 
been largely ignored in the literature. They are, however, crucial if one wants to determine 
whether |$x j i s indeed an honest interaction or only some auxiliary object corresponding 
to the first few terms in an asymptotic series. Our approach is designed to settle such issues. 
In fact, it yields complete mathematical control over the transformed interactions (and the 
unitary transformations giving rise to them): 

(a) We obtain explicit bounds for the error terms, and show that each transformed inter- 
action | i s weu defined in the strong-coupling regime. (In fact, it is exponentially 
summable.) We provide explicit estimations of the radius of convergence, to(n), and 

(n) 

of the exponential decay of the terms <& x in the diameter of X. 
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Figure 1: Zero-temperature phase diagram for H \ = U J2 x e\ n x+ n x - — l J '+J2 x <aA n x+ — 
(J>-J2xeA n x-- Ground states are defined by spin configurations, which at each site, are as 
denoted in curly brackets 
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Figure 2: Shaded region to which we restrict most of the analysis of the generalized Hubbard 
model (P . 
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(b) We provide explicit algorithms to construct the full interaction not only its 

first few orders (restricted to the lowest band). 

As we work in terms of interactions, all the expressions, bounds and estimations are, of 
course, independent of the volume. For the Falicov-Kimball regime, our analysis implies the 
various estimations obtained, by more cumbersome means, in previous studies of the model 



|H [jj, [RJ H S HI (see also the review (TJ). 

Similar calculations, that we leave to the reader, can be done in the region of the plane 
of chemical potentials close to the origin. The low-lying band is now formed by states with 
no doubly occupied sites. This yields a family of transformed interactions } whose 
restrictions to the lower band yield the well known perturbation series often associated with 
the keyword "Gutzwiller projections" . In fact, these series, which start with the very popular 
t-J interaction plus (usually neglected) three-site terms, are the ones usually reported in the 
physics literature (see e.g. f2(J, ||, [7|, |3l|). The half- filled case discussed above is obtained after 
a further projection onto the half-filled band. Our formalism has, in relation to previous 
publications, the advantages (a) and (b) described above. 

(II) Phase diagram at low temperatures 

We are able to infer low-temperature properties of the model ( |1 . 1| ) , in the strong coupling 
regime, in the shaded region of Figure namely away from the manifolds of coexistence of 
ground states, and at and near the segment along the line /i + = as long as we stay clear 
of its infinitely degenerate endpoints. As we will show later by using contour arguments, we 
restrict our attention to dimensions d > 2 and hopping amplitudes t^ that are periodic with 
respect to translations of the sites x, y. Of course, rigorous control of the ground states and 
low-temperature equilibrium states on the "unperturbed" coexistence line u. + = //_ with the 
help of contour expansions will be achieved only for the models without continuous SU(2)- 
symmetry. We shall prove the following results: 

(R3) Inside the (open) regions where H \ has a unique ground state, we apply the cluster 
expansion developed in []8], Section 6] to conclude that the ground state is stable for 
small hopping amplitudes and at low temperatures. That is, in these ranges, one finds 
an equilibrium state for which the observables have expectations close to their values 
for the corresponding ground state. Such a state can be constructed as the limit of 
a sequence of states associated with finite boxes with a ground state configuration as 
boundary condition, and the cluster expansion allows for a visualization (in a very 
precise sense) of the state as "ground state + quantum and thermal fluctuations" . 



(R4) For the Falicov-Kimball regime, \t [ f\ « \e xi \ , we can apply the Pirogov Sinai ap- 
proach developed in || U in conjunction with the transformed interactions j^x"*} de- 
scribed above, to analyze the shaded region of Figure | at or close to the line fi + = 
for strong Coulomb coupling U. In this way we recover and extend a number of pre- 
viously published results, prove some expected features of the low-temperature phase 
diagram and open the way for further systematic analysis of the latter. 

In more detail, the results alluded to in (R4) follows. 




Figure 3: Region (shown shaded) where low-temperature properties of the model ( |1 . 1| ) are 
determined in the strong-coupling regime 
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Falicov-Kimball model without flux (t^ = t for all x, y) 

(R4.1) As soon as t ^ 0, the quantum hopping breaks the infinite-degeneracy present on the 
line < /i + = /!_ < U . Instead, two Neel-ordered ground states appear on this line, 
and in its neighborhood, and remain stable for low temperatures. For the Falicov- 
Kimball model itself (i.e. when t+ x ' is exactly zero), this fact has been proven in p6| 



for chemical potentials on the line itself, and in |30| for potentials in its vicinity. Our 
methods show that the same is true even if one adds some small ionic hopping (a fact 
that also follows from the treatment in [j32]0 or, indeed, any other translation-invariant 
quantum perturbation, which could be long-range but must decay exponentially; see 

(Pi- 



(R4.2) By considering the transformed interaction, |$^|, we can actually show that, for low 
temperatures and t/U small, a phase diagram as in Figure f| appears for the Falicov- 
Kimball model: There are two "strips" of width of order t 4 /U 3 , located around the 
lines /i + — /i_ = ±4t 2 /U for zero temperature, and some smooth deformation [Section 
3.2.1] of them for small nonzero temperatures, where nothing can be said. Between 
them, the states are Neel ordered, while towards the upper left the state is (a thermal 
perturbation of the ground state) all- "+" , and to the bottom right it is all- "— " . Such 



a phase diagram has already been established at zero temperature |To], Pq| , but the 
proof of its stability for small temperatures escaped the methods available up to now. 
Recently, however, this stability has been established by contour methods, related to 
ours, but replacing our partial diagonalization by a study of "restricted ensembles" 



(R4.3) Resorting to the next order in the transformation of interactions, i.e., to |$^|, we 
can see in more detail what goes on inside the two excluded regions discussed above. 
For d = 2, each region splits into four "strips", now of width of order t 6 /U 5 , that 
require further investigation (Figure [5|), and, between them, a finite number of states 
appear, which are quantum and thermal perturbations of certain configurations with 
different proportions of "+" and "— " particles. At zero temperature, these features 



have been formally calculated in [Tj| and rigorously established in |2f| and, by a cluster- 
expansion method, in |32|]. Our formalism shows that they are also present at small 
nonzero temperature, except that the excluded manifolds suffer a smooth (and small) 
deformation. 

(R4.4) At the next order, t 6 /U 5 , each of the eight precedent strips is expected to open up into 
eight finer strips outside of which there is only a finite ground state degeneracy. Our 
formalism would then yield stability of this structure when the temperature is increased 
or further small quantum perturbations are added. In fact, a "devil's staircase" is 
expected to appear in this way, to successively larger orders t n /U n '~ 1 . While we do not 
work out any details we think that this paper provides enough tools for a motivated 
reader to pursue these issues: It provides algorithms for the higher-order transformed 
interactions j, the certainty that the remainders are rigorously controlled to the 
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Figure 4: Phase diagram at zero and low temperatures of the Falicov-Kimball model to 
order t 2 /U, as a function of h — fA+ — /^_. The labels S + , S- and S c b refer to states that 
are quantum and thermal fluctuations of the indicated configurations or its translations 
(open circles correspond to "+" particles and closed circles to "— " particles). Thick lines 
represent excluded regions of width Oit^/U 3 ) located around h± = ±4t 2 /U + 0(ept). The 
small correction ep t is discussed in Sections [4.4| and [5] below 
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Figure 5: Phase diagram at zero and low temperatures of the Falicov-Kimball model to 
order t A /U 3 , as a function of h — //+ — /^_. Besides the states of Figure f|, there appear 
those corresponding to quantum and thermal fluctuations of rotations and tranlations of 
the configurations Si (depicted) and Si (obtained from the Si by a • < — > o interchange). 
The excluded regions (thick lines) have width 0(t 6 /U b ) and their location are determined 
by the values h x = -At 2 /U - 4t 4 /t/ 3 + 0(e#), h 2 = -At 2 /U + 16t 4 /C/ 3 + O(e^), h 3 = 
-At 2 /U + 48t 4 /C/ 3 + 0(ept) and h± = -4t 2 /U + 84t A /U 3 + 0(e#) 
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specified order, and conditions to predict stability. Of course, in order to get precise 
results, one has to come up with an accurate determination of the ground states of 
the classical leading part of the transformed interactions. This tends to be the hardest 
part of the work. The beautiful method introduced in p5| (see also the application in 
HlS| , Appendix A]) could be very helpful at this point. 

An advantage of our methods is their robustness with respect to perturbations: All the 
results sketched above remain valid if we add an arbitrary, but small enough, quantum 
perturbation that respects the symmetries connecting the different ground states (typically 
invariant under translations and 90°-rotations). In particular, they extend to the Falicov- 
Kimball regime (i.e. small, but not necessarily zero), and to models including additional 
perturbations, which can possibly be long-range but must decay exponentially; see ( |2.2(J| ). 

Our methods are not suited for the study of the phase diagram of the S77(2)-symmetric 
Hubbard model, because continuous symmetry breaking is accompanied by the appearance 
of gapless modes. In particular, we cannot apply our phase-diagram technology to the 
effective Hamiltonians that would be obtained via Gutzwiller projections (t-J interaction + 
higher orders). But we emphasize that the perturbation technique (I) is applicable in such 
situations. 



Falicov-Kimball model with flux [t^ = texp(i9 yx ), with 6 yx e (0, 2tv)] 

The previous results can be extended to the model with flux, starting from the study of 
ground states presented in |1S| . For this case, the results (R4.1) and (R4.2) remain valid, 



because the properties involved are independent of the flux. Regarding (R4.3), the positions 
of the excluded regions do depend on the flux, but the ground states around them do not 
(see Section 3.2.1 in Our methods show that these ground states remain stable at small 

temperature and/or under the addition of small translation- and rotation-invariant quantum 
perturbations (e.g. ionic motion). 



1.2.2 3-band Hubbard model 

As an illustration of the use of our perturbation technique for models whose leading in- 
teraction is not strictly on-site, we consider the 3-band Hubbard model, which was first 
introduced by Emery to describe the behaviour of the charge carriers in the Cu02 planes of 
cuprates exhibiting high-T c superconductivity. The model is defined on a two-dimensional 
lattice, a unit cell of which contains one copper and two oxygen atoms (Figure ^j. 

The copper atoms form a square lattice, and there is an oxygen atom between each 
nearest neighbor copper-copper pair. 

The charge carriers are spin-1/2 holes which can hop between the (i-orbital of a copper 
atom and the p-orbital of an adjacent oxygen atom. We shall refer to this as the hopping 
of holes between a copper site and an oxygen site. The Hamiltonian governing the model 
defined on a finite lattice A := A U B, where A denotes the sublattice of copper atoms and 



Tight-binding models 



February 7, 2008 



16 



o o 



o • o • o 



o o 



o • o • o 



o o 



Figure 6: Lattice for the 3-band Hubbard model. The open circles denote oxygen sites while 
the closed circles denote copper sites 

B denotes the sublattice of oxygen atoms, is given as follows: 



H A = \u d J2< + n d x _ + U p J2n p y+ n p y . + U pd £ 

>■ tCA ii f= A A 



<n p v 



xeA j/SA {xy)eA 

x£A yeB xEA ,y€_B 



=: Foa + ^aM. (1.2) 

The number operator for a hole at the site x G A is = (nf + + where nf CT = d\ a d xcn 
for cr G {+,—}, with d xa and rf xo - denoting fermion creation- and annihilation operators for 
a hole of spin a at the site x. The corresponding number operator for a hole at the site 
y G B is = (n^ + + where n v ya = p ya p ya , with and p ya the fermion creation- 
and annihilation operators for a hole of spin a in the p-orbital of an oxygen atom. The 
coefficients U d and U p denote the strengths of the on-site repulsions between two holes of 
opposite spin when they occupy a copper site and an oxygen site, respectively. In addition 
to this on-site interaction, two holes experience a Coulomb repulsion of strength U pd when 
they occupy adjacent copper and oxygen sites. The on-site energies for the copper and the 
oxygen are denoted by the symbols Ed and e p , respectively, and the charge transfer gap is 
given by 

A:=e p -e d . (1.3) 

The hopping amplitude of a hole between a copper site and an oxygen site is denoted by 
t pd and is proportional to an overlap integral of atomic orbitals. Often a small additional 
term representing the direct hopping of a hole between the p-orbitals of two oxygen atoms 
is included in the Hamiltonian, but we neglect it here. 
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Band-structure calculations of cuprates exhibiting high-T c superconductivity show that 
it is not unreasonable to study this model assuming 

A > ; U d > U p > U pd > ; U d > A ; t pd > 0, (1.4) 



and considering the strong coupling regime U d 3> t pd . (For instance, in fl2~2", |23fl, the estimated 
values are A = 3.6eV, U d = 10.5eV, U p = 4eV, U pd = 1.2eV and t pd = 1.3eV.) We observe 
that: (i) we can treat V\(t pd ) as a perturbation of H \; (ii) since A > 0, holes prefer to 
reside on copper sites rather than on oxygen sites, and (iii) since U d is positive and large, 
double occupancy of holes at a copper site is energetically unfavorable. 

We restrict our attention to the situation in which the total number of holes, Nh, in the 
lattice A is equal to the total number of copper sites \A\: 

N h = \A\. (1.5) 

For this choice, a ground-state configuration of if A has a single hole at each copper site, 
there being no holes at the oxygen sites. [Equivalent ly, we could introduce suitable "chemical 
potentials" and work in the region of parameter space where the ground states have this 
property] The hole at a copper site can, however, be either spin-up or spin-down. Hence, 
the ground state of H \ has a 2' A l-fold spin degeneracy. 

By resorting to our perturbation expansion, we obtain an equivalent interaction that only 
involves copper atoms, which includes an antiferromagnetic Heisenberg term. This does not 
settle the question of whether the hopping term V\(t pd ) lifts the spin-degeneracy, because, 
as mentioned above, our methods to determine stable phases do not work in the situation 
where a continuous symmetry [here SU(2)] can be broken spontaneously. Nevertheless, our 
methods show that the transformed interaction is mathematically well defined and provide 
error estimates which should be valuable in future studies. The form of the transformed 
interaction is well known. (Indeed, the fact that the leading terms coincide with those of 
the transformed interaction for the one-band model is an argument favoring the study of the 
latter rather than of the more cumbersome 3-band model). 

It should be noted that the perturbative analysis of models with classical interactions 
that are not on-site is more difficult, technically, than the one of models with single-site 
classical interactions, such as the single-band Hubbard model. Our presentations in Sections 
H and | will serve to illustrate this point. 



2 The Perturbation Technique 
2.1 Foreword 

In quantum mechanics, Rayleigh-Schrodinger perturbation theory can be formulated as the 
perturbative construction of unitary conjugations transforming a Hamiltonian of the form 
H — Hq +tV into a new Hamiltonian, H' = H' Q {t) + V'(t), with the property that H'^t) 
is block-diagonal (in the eigenbasis of H ), and V'(t) is a perturbation of order t n+1 . The 
operators H' Q {t) and V'(t) are given in terms of power series, and it is important to analyze 
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the convergence properties of these series. Kato's analytic perturbation theory p4| , p6[ 
offers criteria that guarantee convergence: Suppose we attempt to block-diagonalize the 
Hamiltonian H with respect to a spectral projection of H corresponding to a bounded 
subset of the spectrum of H separated from its complement by an energy gap 5 > 0. Let us 
assume that the perturbation V is relatively bounded with respect to Hq, e.g. in the sense 
that there exist finite (non-negative) constants a and b such that 

\\Vj>\\ < a||fZoV>||+&|MI, (2.1) 

for all vectors ip in the domain of H . Then perturbation theory for the block-diagonalization 
of H with respect to the given spectral projection converges, provided 



t|a < 1, ^-(a + b)<C, (2.2) 



where C is a constant that depends on the choice of the spectral projection. 

Unfortunately, this result cannot be applied directly to the study of systems of quantum 
statistical mechanics. The reasons are as follows. 



(SMI) First, we do not just want to block-diagonalize a single Hamiltonian H but a family 
of Hamiltonians {-£^a} indexed by an increasing sequence of finite regions A of an 
infinite lattice. In this setting, Kato theory is not applicable, because the norms 
H-HqaII an d ||VaJ| are typically proportional to the volume |A|, and this usually implies 
that the constant b in is proportional to the volume |A|, too. As a consequence, 
the radius of convergence of plain Rayleigh-Schrodinger perturbation series tends to 
zero, as |A| tends to oo. This makes it impossible to establish facts valid in the 
thermodynamic limit by straightforward use of Rayleigh-Schrodinger-Kato analytic 
perturbation theory. (Even if if A has a finite ground state degeneracy, uniformly in 
A, a naive application of Kato theory does not prove that the perturbation series for 
the perturbed ground state energy density, for example, converges uniformly in A.) 

(SM2) Second, to study the phase diagram of the system, the transformed Hamiltonians 
H'(t) must be expressed again in terms of local interactions. That is, the unitary 
conjugations must lead to a block-diagonalization at the level of local interactions, 
rather than just at the level of Hamiltonians. To satisfy this requirement, we find 
that these conjugations must themselves be defined by (exponentials of) sums of local 
operators. Furthermore, the original interaction and the local operators defining the 
conjugation must decrease exponentially in the size of the lattice region in which they 
are localized. 

Keeping in mind these requirements, it is natural to search for transformations whose 
generators S are sums of local terms, S = Sy, where each Sy is localized in the region Y, 
in the sense that it commutes with local operators $x> provided X and Y are disjoint. To 
achieve this, we incorporate two crucial ingredients in our expansion: 
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(11) Defining local ground states (of the leading part H \), we decompose each term <3>x 
according to how it acts on states that are local ground states only in the vicinity of 
X. This is achieved by introducing suitable local projectors (P Y , P Y and P Y below). 

(12) We introduce protection zones to decouple sufficiently far separated excitations. If 
we normalize the H -energj of local ground states to be then the H -energy of a 
configuration of local excitations (local deviations from ground states) separated by 
protection zones is additive in the ifo-energies of the excitations. 

Ingredient (12) is irrelevant for the discussion of models, such as the one-band Hubbard 
model, where the leading interaction is on-site. In such models, the H -energj of excitations 
localized at different sites is additive. But this is untypical, and in the study of more general 
models (12) plays an important role. Regarding ingredient (II), one finds that it has been 
used by some authors || [14], |31], |TTJ, but it has been ignored by others |5], |27|, |39], 0, |(J. 



The latter use "global" projections involving all sites in a given region A, which lead to 
volume-dependent expansions with a shrinking convergence radius as mentioned in (SMI). 
Nevertheless, they appear to obtain correct perturbative results for certain infinite-volume 
expectation values and for (intensive) thermodynamic quantities, like ground-state energy 
densities. The problem with their methods is that they do not offer a handle to proving 
(rigorously) that, indeed, their predictions are correct or to determine whether the corre- 
sponding series are convergent rather than just asymptotic. Expansion schemes such as ours, 
yielding volume-independent effective interactions, are suitable to answer such mathematical 
questions. To accomplish this, they must, however, be supplemented with an appropriate 
construction of ground- and Gibbs states for the transformed Hamiltonian. This is not an 
easy matter, in general, and it is here that quantum Pirogov-Sinai theory will enter the 
scene. 

We are able to provide a complete picture of the models analyzed in Section |^, below, 
where Gibbs states are constructed through quantum Pirogov-Sinai theory. In regions of 
stable phases, the cluster expansion methods we use imply the analyticity of expectations 
and of thermodynamic potentials as functions of the couplings. For those models, therefore, 
we are able to prove that the series obtained for expectations of local observables and for 
intensive quantities are well defined and convergent. 



2.2 Notation, definitions and preliminary results 

In this section, we review the different elements of our general perturbation technique. (For 
further details see ||.) 



2.2.1 Bases and operators 

The Hilbert space. We consider a quantum-mechanical system associated with a finite 
subset A of a //-dimensional lattice ' . To each site x in the lattice is associated a Hilbert 
space TC X . We require that there be a Hilbert space isomorphism (p x : TC X —>■ 7i, for all 
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Figure 7: Spiral order in Z 2 

x G H? , H being a fixed, finite-dimensional Hilbert space. The Hilbert space of the region 
A is given by the ordered tensor product space 

Ha = ®x&aH x . (2.3) 

To avoid ambiguities in the definition of the tensor product (|2.3|) , we choose a total ordering 
(denoted by the symbol -<) of the sites in ~E. V '. For convenience we choose an order with the 
property that, for any finite set X, the set X := {z G Z 1 ', z -< X} of lattice sites which are 
smaller than X, or belong to X, is finite. For example, for v = 2 we can use the spiral order 
depicted in Figure (0). 

Tensor product basis of Ha- Let / be an index set and let {e 3 -}j- e j be an orthonormal 
basis of H. Then {ej}j e /, where e| = ip^ej, is an orthonormal basis of TL X . 

A configuration w on A is an assignment {j x (u)} x( zA of an element j x G / to each x G A. 
For X C A, let uox denote the restriction of the configuration uo to the subset X. The set of 
all configurations in A is denoted by 

C A :={u\j x (u) el,xe A} . (2.4) 

There is a one-to-one correspondence between configurations a; in A and basis vectors 

e{u) := ® x eA(%,( u) (2.5) 

of TCa- A tensor product basis of Ha is given by {e(u;)} a;eCA . 

Local algebras of gauge— invariant operators. Let B(H A ) be the algebra of all bounded 
operators on Ha- Let {U(8) \ 9 G R} be a one-parameter unitary group on Ha (gauge group 
of the first kind) with the property that the vectors e(u) defined in (3.5) are eigenvectors of 
U{6), for all 6. For any X C A, we define a local algebra Ax C B(H A ) of gauge- invariant 
operators with the following properties: 

• If {e(u;)} we c A is an arbitrary tensor product basis of Ha then 

(e(w'),oe(w)) = 0, (2.6) 

unless uj'\a\x — w\a\Xi f° r a ^ operators a G Ax- (Here (ip , cp) denotes the scalar 
product of two vectors ip and ip in Ha-) 
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• U(9) aU(6)* = a, for all a G Ax, for arbitrary X C A; (gauge-invariance). 

• If X C Y then ^ C Ay. 

• If X n Y = then 

[a,6]=0, (2.7) 

for any a G Ax and any 6 G Ay. 

Example: Gauge-invariant polynomials in fermion creation- and annihilation operators. Let 

{c xa , 4 CT , | x G A, a = +, -} (2.8) 

denote the usual fermion creation- and annihilation operators satisfying the canonical anti- 
commutation relations 

{<W, c ya ,} = {c4, c\ a ,} = 0, (2.9) 

and 

{Cxa, <y ya >} = $xy <5aa>- (2.10) 

Of course, a is the spin index. We define 



U{0) = exrj(i9N) , where N = ^ c\ a c X(J . (2.11) 
The local algebras Ax, X C A, are defined by 



16A 
^6{ + ,-} 



Ax = {a | a is an arbitrary polynomial in c xa , c^, x G X, a E {+,—}, 

with the property that U(9) a U(6)* = a}. (2.12) 

The algebra of local, gauge- invariant operators. Since Ax C Ay, for X C Y, we 

may consider the (inductive) limit Vx /r -^x- The algebra A is defined to be the closure 
of Vx /-z? Ax in the operator norm and is called the algebra of all local, gauge-invariant 
operators. 

2.2.2 Interactions and ground states 

Classical operators and -interactions. An operator a is "classical" w.r.t. a given tensor 
product basis {e(a;)} if and only if e(u) is an eigenvector of a, for all uj G Ca- An interaction 
{$ox} is said to be a "classical interaction" (w.r.t. {e(u;)} a;gCA ) if and only if the following 
conditions hold: 

1. $ox ^ <Ax is a classical operator w.r.t {e(u)}, for all X C A. 

2. More precisely, 

<$>oxe(uj) = $ x(w)e(w) , $ x(w) G R, (2.13) 
where $ox(^) only depends on u\x- 
[Note that 2. follows from 1. and from fl2.7| )l. 
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m-potentials. Let $o = {^ox} be a classical interaction. We shall always assume that 
there is at least one configuration, ujq, minimizing all $ox ; i.e., 

%x (vo) = min $ ox (u), (2. 14) 

for all X; (this assumption holds if $ is given in terms of a so-called m-potential pi|). We 
normalize the operators <3>ox such that $ x(^o) =0, for all X. Thus 

$oxM > and min$ x(^) = 0, (2.15) 

for all X. 



Local ground states. A configuration u is said to be a local ground-state configuration 
for a classical interaction <3> an d a subset X of the lattice if 

$oyM = for all Y C X. (2.16) 

The local Hamiltonians H x — J2ycx ^oy thus have the property that 

H ox e(u) = 0, (2.17) 

whenever to is a local ground state configuration for X. A state e is said to be a local ground 
state for a region X if e is an arbitrary linear combination of the vectors e(u), where uj 
ranges over the local ground state configurations for X. 

We shall always assume that the spectrum of H x has a gap above which is bounded 
from below by a positive X-independent constant, for all bounded sets X. This assumption 
can be derived from suitable assumptions on the classical interaction $ (^o should be an 
m-potential with a finite number of periodic ground states, see ||21|| ). 



The full interaction. We consider quantum lattice systems whose dynamics is encoded 
in an interaction $ which is assumed to be of the form 



$ = $o + Q(A), (2.18) 

where $o = {^ox} is a finite-range classical interaction given in terms of an m-potential, 
and Q(X) = {Qx(A)| is a perturbation interaction with A = {Ai, . . . , A^} a family of real or 
complex perturbation parameters. For brevity, we shall say that a set A is a classical bond 
if $ x 7^ 0, and a quantum bond if QxW ^ 0. The operators $ x an d Qx(A) belong to the 
local algebra Ax- A polynomial p(X) is said to be of degree n if 

p(\) = Y, C n,..n k K 1 ---K k , (2-19) 



with c ni ...„ fc = whenever J2j=i n j > n - We assume that the leading order in A of the 
interaction Q(X) in (|2.18|) is of degree 1 (i.e., linear in A). 
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We also assume that the interaction Q(X) = {Qx(X)} is either of finite range or decays 
exponentially in the size of its support, i.e., that there exist an r > such that 

E vMQx(A)Qx(A)K s(x) < oo, (2.20) 

where s(X) denotes the cardinality of the smallest connected subset of the lattice which 
contains X. The interaction Q(X) will be treated as a perturbation. 

Remark. Roughly speaking, we shall always split the Hamiltonian H\ into a sum of 
an unperturbed operator H Q \ and a perturbation Va(A) in such a way that all low-lying 
eigenstates of if a, corresponding to eigenvalues separated from the rest of the spectrum 
of if a by a positive (A-independent) gap, correspond to degenerate ground states of iioA- 
More precisely, we choose the classical interaction $ i n such a way that all local low-energy 
eigenstates of the local Hamiltonians H x, for X C A are exactly degenerate in energy; 
small degeneracy-lifting terms are systematically put into the interaction Q(X) defining the 
perturbation V\(X). In a general exposition of our methods, this is a very convenient way 
of ensuring that ifo-energies of sufficiently far separated local excitations are additive. (Of 
course, in a concrete model, additivity of ifo-energies of local excitations may hold for 
independent reasons; see e.g. Section ||) 

2.2.3 Projection operators. Protection zones 

"Diagonal" and "off-diagonal" operators. Given a partition of unity 

1 = £P, , P i P j = 5 tj P j , P; = P,, (2.21) 

i=i 

and an operator Q, we set 

Q ij I'.QI'r (2.22) 

We will call the operators Q u "diagonal", and the operators , i ^ j, "off-diagonal" (with 
respect to the given partition of unity). 

Protection zones To deal properly with the local nature of different operators we define, 
for any x G Z^,a so-called R-plaquette centered at x : 

W x := {y e TL V : \ Vi - Xi \ < R, fori < i < u}, (2.23) 

where R is the range of the interaction $o- For any finite set X C 7L V ', its covering by 
P-plaquettes is denoted by 

Bx := U xeX W x . (2.24) 
The set Bx \ X can be interpreted as a protection zone around X. 
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Local projections. We introduce some special projection operators on the Hilbert space 
Ha: 



1. If Y C A then P Y is the orthogonal projection onto local ground states for Y. If Z C Y 
then P| D P£. 



2. The orthogonal projection onto the space of states which are ground states on Bx \AT, 
but fail to be ground states on X: 

Pb x '■= Pb x \x ~ Pb x - (2.25) 

Hence Pi projects onto the space of states which have an excitation localized in X. 
The set B x \ X acts as a "protection zone" introduced to ensure the additivity of 
energies of disconnected excitations [see fl2.30|) below]. 

3. The projection onto states with excitations in the "protection zone" Bx \ X: 

P 2 b x ■= 1 " P°b x \x, (2-26) 

where 1 is the identity operator. 

Additivity of excitation energies. For each finite X C Z^, we can decompose the 
Hamiltonian H \ as follows: 

Hoa = H ox + Hqb c x + Hob x \x, (2.27) 

where B c x = A\B X , 

Hoa = E $ or- (2.28) 

YcA 

and 

H 0A = E $ oy (2-29) 

If ip G Ran Pb x \ x tnen 

HoaV' = (H x + Hob x \x + H 0B c x )ip 

= (H x + H 0B o x )t/j. (2.30) 

This follows, because HoB x \xip = H b x \x PbxXx^ = 0- Hence, on states ip corresponding 
to ground-state configurations in the region B x \ X, energies are additive. 
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2.2.4 Operator identities 

The Lie-Schwinger series With this expression we shall refer to the following expansion: 

1 



e A Be- A = B + [A,B] + -[A,[A,B]] + ... 

OO 1 

= £-TacfM(5). (2.31) 

n=0 U - 

Here we use the notation 

adA(B) = [A, B] ; ad 2 A(B) = [A, [A, B]] ; ad n A(B) = [A, ad" 1 ' 1 (B)} (2.32) 
and the convention 

ad° A(B) = B . (2.33) 

The "ad -1 " operation Given a self-adjoint operator H whose spectrum consists only of 
eigenvalues Ei, % — 1,2, . . ., and an operator Q which is purely "off-diagonal" with respect 
to the partition of unity given by the projections P % onto the eigenspaces corresponding to 
each Ei, we define 

ad-^QH^TP*— 9— (2.34) 

The right hand side is well defined, because Q is an off-diagonal operator. Of course, ad -1 
is the operation inverse to ad, i.e., 

[H, ad' 1 H(Q)] = Q. (2.35) 
2.3 First-order perturbation theory 

The first step in our perturbation technique consists in eliminating "off-diagonal" operators 
of lowest order in A. Here "off-diagonal" refers to the partitions of unity 

l = Pl x +P l Bx +Pl x - (2.36) 

In the following, we usually suppress the explicit dependence of the operators Qx(^) on A. 
We rewrite Qx = Qx(A) as 

Qx = QZ+Qb x +Qb x , (2.37) 



^°3 X + Q°B X + Qb x : 

)00 :„ „ u r 



where Qb x is a "diagonal" operator defined by 

Qb x ■■= Pb x QxP° Bx , (2.38) 
Q B 1 X is an "off-diagonal" operator given by 

Qb\ ■= Pb x QxPL + PLQxPL, (2.39) 



Tight-binding models 



February 7, 2008 



26 



and Qb x is the remainder 

Qk ■= Pl x QxP 2 B x + P l B x QxPi x . (2.40) 

Note that, since Qx £ Ax, the operators P 1 b x QxPb x an d Pb x QxP % b x vamsn , for i = 0,1. 



Using the decomposition ( 2.37|) we can write the Hamiltonian if a as follows: 

H A = H 0A + V A 00 + V* 1 + V R , (2.41) 

where ^ A 00 = £ XcA Qb x , V^ 1 = ExcA Qb\ , and V? = ExcA Qb x ■ 

There is a slight subtlety, at this point, related to boundary conditions: The definition 
of the operators P Bx , i = 0, 1,2, and of Q°b x , Qb\i Qb x ^ oy regions X with the property 
that E>x H A c is non-empty must, in general, be modified in such a way that the boundary 
conditions imposed on the configurations on A c are properly taken into account. We shall 
not enter into a detailed discussion of this (primarily technical and usually straightforward) 
issue. In fact, we shall usually think of periodic boundary conditions for which the issue 
does not arise. 

For notational simplicity we henceforth suppress the subscript A. Following the guidelines 
explained in the Introduction, we search to eliminate the first-order off-diagonal terms QbK > 
through a unitary transformation C/W(A) = exp(Si(A)), where Si (A) is of degree 1 and is 
given by a sum of local operators: 

S 1 (\):=J2S 1Bx (A). (2.42) 
x 

In the sequel we shall also suppress the explicit A-dependence of the operators Si (A) and 
Sib x (A). In order to gain mathematical control of the resulting expressions, it is essential 
that each Sib x be a local operator. 

By the Lie-Schwinger series (|2.31|) , the unitary operator yields the transformed 
Hamiltonian 

#M : = e Sl He~ Sl 

= H + V 00 + V 01 + V R 

+ Y -4ad n S a (H + V 00 + V 01 + V R ). (2.43) 

To leading order, the operators Si, V 00 , V 01 and V R depend linearly on the perturbation 
parameter A. We wish to eliminate V^ 01 . This leads to the condition 

ad# (Si) = V 01 . (2.44) 

If this equation is satisfied, the Hamiltonian (|2.43|) becomes — singling out the terms up to 
second order in A — 

#W = H + V 00 + V R + V 2 

+ + V 01 ), (2.45) 

^ 2 n\ Ti + 1 
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with 

yOl 

V 2 :=zdS l (V 00 + V R + —). (2.46) 

Condition ( |2.44| ) reads 

Y,*dH (S 1Bx ) = J2Q° B l x , (2.47) 
x x 



which, given (|2.35|) , suggests the choice Sib x = ad 1 Hq (Q°b x )- The locality of the operator 
S\b x follows from ( |2.30|) , which implies that 

ad- 1 fl (Qg f ) = ad-^ox (Q £ x ) ■ (2.48) 

This is because the operator Qg 1 vanishes on all states which are not local ground states 
for H b x \x (see ( |2.6p and ( |2.17| )), whereas H 0Bx measures the energy of the configuration 
outside X; see definition fl2.29|) . Hence, (|2.47|) is satisfied if we choose 

S 1Bx = ad-^ox (Q°£ x ) . (2.49) 



Selfadjointness of H ox and Q B l x implies anti-selfadjointness of S±b x and hence of S±. The 
identities ( p. 42 ) and ( P-49| ) show that Si is given by a sum of local operators. The condition 
that Sib x should be a local operator (i.e., Si Bx G Ab x ) motivates our introduction of the 
"protection zone" B x \ X. In || it is shown that the family {Sib x } satisfies the required 
summability condition if the initial interaction {Qx} does. 

More precisely, if P % denotes the projection onto the eigenspace of H ox corresponding to 
an eigenvalue E { , E = < E x < E 2 < . . ., then from (|2.34j) 

noi 

For Qb\ as i n ( P-39| ), i ^ j (and either i = 0, or j = ); thus E{ — Ej ^ 0. Hence S\b x is a 
well defined operator. Let us point out that, with these definitions, 

00 , n>R i ^%2\ 



v 2 = j: ^s 1Bxi (q°b X2 + Qb X2 + 



x l,X 2 

B Xl nx 2 ^ii 



E V 2BxiUX2 . (2.51) 



:S x 1 ux 2 ■ 
x 1 ,x 2 
B x 1 nx 2 ^di 



The transformed Hamiltonian ( |2.43|) can be written in the form 

H (i) = K (i) +jR (i) ; (2.52) 

where the operator K^ 1 ' contains the leading-order block-diagonal contributions: 

=H + V 00 + V R + V 2 00 + V 2 R , (2.53) 
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where the last two terms are sums of operators 

:= PR V 2Bx , 1V PS (2.54) 

^B Xl ux 2 &x 1 ux 2 za x 1 ux 2 B Xl ux 2 v ' 

V 2B Xl ux 2 := P B XiUX2 V 2B XiUX2 Pb XiU x 2 + Pb XiUX2 V 2B XiUX2 Pb XiUX2 - (2.55) 

This leading part K^ 1 ' basically corresponds to the first-order expressions usually presented 
in the literature on perturbation expansions. The "remainder" RP-\ on the other hand, 
starts with the leading "non-diagonal" terms and includes all higher-order contributions 
arising from the Lie-Schwinger series ( |2.43| ) 

P« = Vf 

+ E E ^ ^S lBxn (adfi^ (... (Qb Xq + Q R Bxo + ^~ Q B \ ) ■ 

n>2 X ,X lt ...X n : "'■ ~ 
X U...UX„= c.s. 

(2.56) 

where X U . . . U X n = c.s. denotes the condition 

Xt n By._, ^ for 1 < i < n, (2.57) 
with Fi_i =X U...U Xi-i, and 

^2B XlUXa := P Bx 1 uX 2 V2B X 1 uX 2 P B XlUX2 + P B XlUX2 V 2B XlUX2 P BxiUX2 - (2.58) 

As shown in ||, the series (|2.56| ) converges, in fact absolutely and exponentially fast, if the 
original interaction Q satisfies ( |2.20| ). 

In fact, the transformed Hamiltonian comes from an interaction 

$ (i) = Q« (2.59) 

where \&g an d are the interactions whose terms yield and P^ 1 -* respectively. Ex- 
plicitly, \E f o 1 ' ) = {^ofiy} where the ^ob y are nonzero only for Y = X or Y = X\ U X2 with 
Bx 1 fl X 2 7^ and X, Xl, X 2 (original) bonds and take the values 



*?i x = ^ x ^p^ + pi x $oxP^ + g°B x + g| x + v 2 ° B o x + ^ x (2.60) 

= +^ . (2.61) 

u£>x 1 ux 2 ^- t >x 1 ux 2 '-"XiL^ v ' 

All the operators involved are "diagonal" with respect to the partition of unity fl2.36j ) or the 
analogous partition with X replaced by X\ U X 2 . The interaction has only terms of the 
form {Qb^}, where Y is a "c.s. "-connected set of quantum bonds: 

Y = U] =0 Xj = c.s. , n > 1 , (2.62) 
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and 



Von n = 1 

~m 2 J B XoUXl 

Qb y = 1 n N (2-63) 

The analysis of || shows that the operators Q% , which make up the remainder i?^ 1 ), 
decay exponentially in the size of their supports. More precisely, there exists a positive 
constant r\ > such that 

E V / tr[gg ) *(A)gg ) (A)]e- s ^ = 0(A 2 ), (2.64) 

Z30 

where s(Z) denotes the cardinality of the smallest connected subset of the lattice contain- 
ing Z. Furthermore, the part of fl2.64Q involving "00" -components is of order 3 or larger 
in A. It is, therefore, natural to expect that this interaction does not contribute to 

the ground-state energy of the transformed Hamiltonian to second order in the per- 

turbation parameter A. This fact is easy to see in finite volume, but, as mentioned before, 
its verification in the thermodynamic limit — at the level of energy densities — requires a 
suitable construction of (infinite volume) ground states. For this reason we shall be able to 
confirm this expectation only in those cases were the transformed interaction satisfies the 
hypotheses of the quantum Pirogov-Sinai theory of Section |4.2| . (This excludes models with 
a continuous symmetry.) 

The methods discussed in Section [|for the study of the phase diagram rely on the choice 
of a suitable tensor product basis {e(u>)} such that the operators ^oB y can be decomposed 
into 

*8t = (2-65) 

where $oB y i s exactly diagonal in the tensor product basis {e(u)} — i.e., $oB y i s "classical" 
with respect to {e(oo)} — and is a perturbation which is usually non-diagonal. Such 
a decomposition is always possible. The optimal choice of {e(u)} is the one that renders 
the perturbation as small as possible. It is worth noting that the new unperturbed 

(classical) interaction is of finite range. 
We finally set 

$W = $W + q(i) j (2.66) 

where 

^ = {$£U (2-67) 

and 

Q W = {Q { Z + n$ Y }. (2.68) 

In the simplest examples, the only contribution to the ground state energy of the transformed 
interaction to second order in the perturbation parameter A, arises from the classical 
interaction ^q 1 '' . In such cases one may gain heuristic insight into the structure of ground- and 
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low-temperature states of $ from the ground states of $q . For example, large degeneracies 
in the spectrum of the Hamiltonian determined by $o ma y turn out to be lifted in the 
spectrum of the one determined by $g • Hence this last interaction, though also classical, 
ought to lead to more accurate insight into low-temperature properties of the system than 
the original $o- The examples below illustrate instances where these expectations can be 
rigorously confirmed. 



2.4 Second-order perturbation theory 

Next, we search for a unitary transformation which removes the off-diagonal part of the 
original perturbation, {Q°b y }y, to order |A| 2 . A possible approach (certainly not the only 
one, see the discussion in Section |2J)| below), is to transform the interaction $ < - 1 - ) so as to 
eliminate the lowest order off-diagonal part of the perturbation This leads one to 

consider a unitary transformation of the form 



where S\ = Si (A) = J2xcA Sib x is defined by ( [2.49D , while S 2 (A) is determined from the 



above requirement. The generator S 2 (A) must have leading terms of degree 2 and is required 
to be a sum of local operators: 

S 2 = S 2 (A) = ]T S 2By (A). (2.70) 

YcA 



We notice that the choice (|2.69| ) does not amount to an iteration of the first-order pro- 



cedure, because we keep using the partitions of unity determined by the spectral projections 
of the original classical part $o ; rather than those corresponding to the transformed $q . 



We comment in Section |2Jj, below, how an iterated method (Method 3) would proceed. 

To alleviate the forthcoming formulas let us suppress the explicit A-dependence of the 
operators Si (A) and S 2 (A). The unitary transformation yields the Hamiltonian 

#( 2 ) := e S2 e Sl H e~ Sl e~ S2 



-s 2 



(2.71) 



where is the Hamiltonian given by (|2.45|) . Applying the Lie-Schwinger series we see 
that the term V 2 01 can be removed by defining 



So 



zd^HoiV, 



oi-. 



= ad 1 ^ox 1 ux 2 (^ / 2°i XlUX2 ) 

x 1 ,x 2 



2Bx 1 ux 2 



x 1 ,x 2 
Xin% #0 



^2.72) 
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With this choice, we obtain — singling out the terms of up to fourth order in A — 

#( 2 ) = H + V 00 + V R + V$° + V 2 R + V 3 + V 4 + T (2.73) 

where 

t/OO yR 2V 01 V 01 

V 3 = ad 2 ^(^- + V — + flj-) + adS 2 (V™ + V* + \), (2.74) 

T/00 yR ot/01 

^ = ad3 ^(^r + ^r + — )' ( 2 - 75 ) 

and T consists of all terms of order n with n > 5. 

The transformed Hamiltonian if ( 2 ) can be written as 

H® = K®+R®, (2.76) 

where 

K <2) = K (l) + yOO + y * + ^00 + yR (2 _ 77) 

is given entirely in terms of an interaction = {^ Q b y }, with the property that the 

(2) 

operators ^ b y are "diagonal" with respect to the partitions of unity 

1 = P°b y + Pb y +P 2 b y - (2-78) 

and are of order A; in A, with < k < 4. All terms which are "off-diagonal" — starting with 
V® 1 + V® 1 — , and all "diagonal" terms of higher orders, are included in the remainder , 
which is given in terms of an interaction = {Qb y }- The analysis of || shows that these 
operators satisfy 

E MQf\W { z\X)Y 2iZ) = 0(A 3 ), (2.79) 

for some r 2 > 0, and that a similar sum involving "00" -components is of order 5 or larger in 
A. 

Further, as in ([2.651) , we can write 



tf( 2 > = + (2.80) 

(2) (2) 

where $q = {$o# y } is a finite range classical interaction, i.e., diagonal in a tensor product 
basis {e(u>)}, and is a generally non-diagonal perturbation. One expects that only the 
classical interaction $q contributes to the ground-state energy to order 4 in A. 



2.5 Synopsis of perturbation theory 

The above procedure can be iterated to any finite order in A. To n-th order, we aim at 
constructing a unitary operator, C/jy(A), such that 

H£\\) = u£\\)H A (\)U$>\xf (2.81) 
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has the property that 

ft = E*i n) (A), (2.82) 

ZcA 

with 

P$, s *fWPh, + Pk^\X)Pk = 0(\M n+1 ) (2-83) 
for all Z. We need to ensure that the conjugation of an interaction with exponential decay 
by the operator uj^\\) is again an interaction with exponential decay. In || we present two 
possible structures for an operator (A) compatible with this requirement (we drop the 
subscript A): 



Method 1: 



Method 2: 



[t ) (A) = n^ (A) ) (2.84) 



f/f(A)=exp£>f (A)), (2.85) 

3=1 



In both cases, 



5j o) (A) = X)5g ) (A) ) « = lor2. (2.86) 

z 

,(a) 



where the local operators S- z (A) are of degree j in A. The idea is to determine these 
operators recursively in such a way that the condition (|2.83 ) holds, for n = 1, 2, 3, . . .. Given 



S^y , . . . , S%y , for arbitrary F, (|2.83|) uniquely fixes the "off-diagonal" contribution to S%9 1Z , 
for any Z and for any k < n. Explicit formulas may be found in M. The method described 



in Section 2.4 follows method 1. 



There is a third way of organizing the recursive perturbative construction of unitary op- 
erators U( n >(X) approximately block-diagonalizing the Hamiltonian H, inspired by Newton's 
method and similar, in spirit, to the renormalization group strategy. Recall that first order 
perturbation theory [see (|2.66|) in Section |2.3|| yields a unitary operator £A 1 '(A) such that 



UM(\)HUW(\)* =: = + V®, (2.87) 

where Hq = J2y ^ob y > = {®()b x j * s a classical interaction, and the interaction giving 
rise to := J2y QbI nas "off-diagonal" terms of order 2 in A. The third method to 
construct the conjugations is based on keeping track of the low-energy spectrum of the 
Hamiltonians determined by {•I'oBy} (rather than by {$ B y }), and, in particular, of their 
ground states. This gives rise to a partition of unity 

1= ^P° Bx + ^Pl x + (1) P 2 b x , (2-88) 

which serves to decompose the operators Qx into 

Q ( x = [QiY + [QiY + [Q { xY, (2-89) 
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as in ( |2.37|) , and hence 



(2.90) 



in analogy to ( |2.41| ). We now proceed as in ( |2.42| ) through fl2.50| ) of Section |2.3| . This yields 
an operator 

Si 3 \\) = J2si 3 l(\) (2.91) 



such that 
has the form 



x 



(2.92) 



(2.93) 



where the perturbation does not have any "off-diagonal" terms of order < 2, with 
respect to the partition of unity (|2.88j). We may now proceed recursively, in the same 



manner, eliminate off-diagonal terms to ever higher order in A. 

The analysis of || remains valid for this case too, and provides proofs of convergence 
and summability of the different interactions. 



3 Perturbation expansion for the one-band Hubbard 
model near half-filling 

3.1 Preliminary remarks 

In this section we derive the effective Hamiltonians for the one-band Hubbard model, defined 



by the Hamiltonian (|1.1|) , near half-filling by applying the perturbative unitary conjugation, 
Method 1, described in Section |2~5| . This model has a number of simplifying features that 
makes it ideal for a non-trivial exemplary application of our methods. The simplifications 
stem from the fact that its leading classical interaction [formula ( [3. 3D below] is a sum of 
on-site terms. First, this implies that every eigenvalue of Hqx,X C A [defined as in ( |2.28| )1, 
is a sum of on-site energies e x (s) [where s denotes the on-site configuration and x G A] which 
for this model are 

e<n(T) = ~H+ , = -fi- , e a (0) = , e x .(t, |) = U - /i+ - n- . (3.1) 

Hence its zero-temperature phase diagram is easy to determine (Figure [I]). 

Second, the classical interaction has range R = 0, and there is no need for protection 
zones: 

B Y = Y, (3.2) 

for arbitrary FcA. 

Below we shall consider, for the whole shaded region of Figure @, projections onto the 
entire band of single-occupancy states. These states are all simultaneously (local) ground 
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states of the classical interaction along the line n + = (in the complement of this line this 
degeneracy tends to be lifted). 

We present the block-diagonalization of the Hubbard Hamiltonian only to second (Section 



3]4j) and third order (Section |3.5|) in the hopping amplitudes. But we believe that the 
following discussion, together with Section 4 of 0, provide enough details for a motivated 
reader to pursue the procedure to an arbitrary order. 

3.2 The original interaction 

As mentioned in the Introduction, it is necessary to implement our perturbation scheme at 
the level of interactions instead of Hamiltonians. The dominant term, Hq\, of the Hamilto- 
nian ( |1.1|) can be expressed in terms of an on-site interaction $ = {$ox } where 

h k 
®ox = Un x+ n x _ - -(n x+ - n x J) - -{n x+ + n x J), (3.3) 

with x G A and 

h := fi + - /c:=/i + + yU_. (3.4) 

As in ( p.28| ), this interaction defines local Hamiltonians 

H 0Y = E • (3.5) 

The quantum perturbation is the kinetic-energy term, and the perturbation parameters are 
the hopping amplitudes t = {t b + ,t b _,t b + ,t b _}, where b = [xy] denotes an ordered pair of 
nearest neighbor sites and b is the same nearest neighbor pair, but with reverse ordering. 
The perturbation is assumed to be translation invariant, that is t±b only depends on the 
direction, y — x, of b. We shall write 

t = max{t b + ,t b _,t b + ,t b _}. (3.6) 
We write the quantum interaction as 

Qx = Qx(t) = Q [xy] (t) + Q [yx] (t), (3.7) 
where X denotes the unordered pair (xy) and 

Q b = Qb{t) = £ Q ab {t) (3.8) 

cr=± 

with 

Qcrb = Qabit) := tc\ a C ya . (3.9) 

Moreover, if the Hamiltonian is assumed to be self-adjoint then 

t\ = (4)*. (3.10) 
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3.3 The local projections 

We start by introducing partitions of unity giving rise to appropriate decompositions of the 
hopping terms in the Hubbard Hamiltonian into "diagonal" and "off-diagonal" operators. 
The most interesting part of the shaded region of Figure |] is the vicinity of the line of infinite 
degeneracy. On this line, the ground states of $o axe (linear combinations of) configurations 
with precisely one (up-spin or down-spin) electron occupying each site. For Y C A, the 
orthogonal projections onto the subspace spanned by the ground states of H 0Y are of the 
form 

P Y := I] P* (3-H) 

with 

P Q x :={n x+ -n x _f. (3.12) 



As remarked above, as the classical interaction ( |3.3|) is on-site, i.e., of zero range, there 
is no need for protection zones. This implies, first, that we can define the projection Pq Y by 
setting 

l = Py+Py, (3.13) 

and, second, that we can use these projections P Y and P Y throughout the whole shaded region 
of Figure 0. The operators P Y are projections onto states with single occupancy. This family 
of states includes, therefore, ground states and low- lying excitations The operators P Y , on 
the other hand, project onto the complementary subspace corresponding to of high-energy 
excitations, formed by states with at least one site in Y occupied by zero or two electrons. 

3.4 First-order perturbation theory 

Let us perform the first step of the perturbative analysis of Section ^|, namely the block- 
diagonalization of the interaction to second order in the perturbation parameters. In this 
section, X is always used to denote a pair of nearest neighbor sites. 

We consider the perturbation defined in (|3.7|) - (|3.9|) and decompose it into the form (|2.37|) 



Q x = goo + goi + qR^ using the pro j ections po anc i pi defined in flSTTTHCT) (P Y = 
for all Y C A). It is simple to check that 

= 0. (3.14) 
Our goal is to construct a unitary transformation, 

= e Sl W, (3.15) 
with the property that H^\t) := U^(t)H/JJ^(t)* only contains off-diagonal operators of 



order 2 and higher in t. Following the general ideas discussed in Section |2.3| , we attempt to 
construct an operator S\(t) of the form 



Si(t) = ££Vx(t) = E E&sw 

X a=+- b 



(3.16) 
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The formulas for the operators Si a b can be inferred from Section |2| [formula ( [2.49 )1: 

S\ a (xy) = ad -1 Ho^y^cl^ c yo ) . (3-17) 

Conjugating the Hamiltonian H of the asymmetric single-band Hubbard model by the opera- 
tor e Sl ^\ yields an effective Hamiltonian [see ( [2.45| )1 given in terms of a new interaction, 
$«, determined by (P32|)-(P39l) (with B x = X). 



Further, as explained in Section ^73] [see (|2.65|) ff.], to be able to apply the Pirogov-Sinai 
theory of quantum phase diagrams, we need to decompose the operators $^ as in (|2.66| ), 
namely in the form 



*$ = *$ + Q ( i\ (3-18) 

(i) 
x 



where {$ x} i s a finite range classical interaction and {Q^} is a perturbation with [(J*- 1 ' 1 ' 1 " 
0(t 2 ). For the Hubbard model, 



^ = = (3,9) 

( \P x aAS lx {Qx)Px if Xis a pair of n.n. 

It is not hard to see that the two-bond terms V x ° lUX2 = ^[&dSiB Xl (Q°b X2 )+Sib X2 (Qb Xi )] 00 are 
zero if X\ ^ X2. The second expression on the R.S. of (|3.19|) involves operators corresponding 
to an electron of a particular spin hopping from a singly occupied site to a nearest neighbor 
site occupied by an electron of opposite spin, thus momentarily creating a doubly occupied 
site and a hole. They then restore the condition of single occupancy of the sites by causing 
one of the electrons on the doubly occupied site to hop to the neighboring hole. 

In fact, the interaction ( 3.1 9| ) coincides with the first-order calculations often presented 



in the literature since the initial observation of Anderson ]]J. It is conveniently expressed in 
terms of the spin operators 

S x := cLSss'Cxs' (3.20) 

s,s'=+,- 

where S = ^(cr 1 , a 2 , a 3 ) and a 1 , a 2 , a 3 are the standard Pauli matrices. Note that 

Sl= l -(n x+ -n x _), (3.21) 

and 

St:=S l x ±iS 2 x = cl ± c XT . (3.22) 

The operators S x , S~ = (S x )* and S 3 form a basis of a representation of the Lie algebra of 
577(2) on the four- dimensional Hilbert space H x = span{0, f, |, f!}. These operators anni- 
hilate the subspace spanned by {0, ||} and act irreducibly (in the spin-1/2 representation) 
on the subspace spanned by {|, !}. In particular, 



pO q3 pO c3 pO c± pO c± (o no\ 

r x ^x r x — D x 1 r x °x r x — °x ■ 

The operators adS^^Q^), with a = +, — , cause the same electron to hop in both steps 
of the process, whereas the operators &dS la .i(Q®} b ), with a 7^ a', cause electrons of opposite 
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spins to hop in the two steps. These processes can be expressed in terms of a sequence of 
creation and annihilation operators of the up-spin and down-spin electrons 

1 



P°x adS^(Q^) P° x 



u Fx 



pO 

r x 



f b f b 

L a L a' P 



u 



C,,„C xcr C iC yrT ' -f- C C vtT C i c xo 



C ya' C xcr' c xa c y<T + C xa' C y^' C ya C xcr 



pO 

r x- 



(3.24) 



For a = a', (^2|) yields 

E P x ^S laX (Q°Jx)Px 



2(KI 2 + I^| 2 ) D o 



U 



P x [n y+ (1 - n x+ ) + (1 - n y+ )n x+ ]P x 



U v ' ■" 4 

To arrive at (|3.25|) , we have made use of the identities 

P x n ya (l - n xa )P° x = P° x (1 - n y{ ^)n x ^ a) P x , 

and 

P X \p>x-\- n y- + n x- n y+] P X 

for X = {xy}. In turn, this last identity follows from ( |3.21| ), ( |3.23| ) and the relations 



pO 

r\ q3 q3 I X 
-ZD x D y i — 



P°rwP° = P°(l-n x( _ a) )P° 



and 



Px i n x+ + ^a;-) Px — Px ■ 

Similarly, for a ^ a', ( |3.24| ) yields 

2 (t\ t b _ + t b _ t b 



(3.25) 

(3.26) 
(3.27) 

(3.28) 
(3.29) 



E & dS lab(Q 



ou 

ah) 



U 



s^s x + s;s: 



4 (t b + t b _ + t b _ t b 

u 



-p 



x 



ql ql _j_ q2 q2 

x y * x y 



A 



(3.30) 



^From ( |3.19| ), ( |3.25| ) and ( |3.30| ) we conclude that $q^- is the well known Heisenberg 
interaction 



(i) 
ox 



-h S x 2 , 



2 fe| 2 +|t 6 _| 2 



pO 



2{t\ t b _+t b _ t\) 
U 



Ql Ql I Q2 q2 
°x D y i" 



if X = {x} 

if X is a pair of n.n. 

(3.31) 
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[h and k being denned in (glM , and = 0, for all other X C Z^. 

If the hopping amplitudes depend neither on the direction of the hopping nor on the 
spin orientation — t±b = t ± i = t — then the original Hamiltonian ( |1 . 1|) is the standard 
one-band Hubbard Hamiltonian, and the interaction ( |3.31[ ) is that of an isotropic Heisenberg 
model. This model is difficult to treat rigorously because of its SU(2) symmetry, and it is 
not surprising that the symmetric Hubbard model lies outside the scope of the methods 

the effective interaction (|3.31 



we discuss in Section [4[ For high asymmetry, 
corresponds to an antiferromagnetic Ising interaction, perturbed by a small spin-flip term 
This interaction exhibits a first-order phase transition to a Neel phase. This applies in 
particular to the Falicov Kimball model, for which t b + = and t b _ = t. 

^From the discussion of Section |2.3| we know that the other terms of the transformed 
interaction are small. The terms corresponding to the rest of called in 02.651) , 
involve only transitions within the excited band, while the "off-diagonal" remainder is of 
second order in t [equation ( |2.64|) 1. 



3.5 Second-order perturbation theory 



Let us now pursue the block diagonalization process to one further order. The objective is 
to find the classical interaction of (|2.80 ). General theory tells us that the remainder is 



small [bound ( [2.791) 1. in analogy to the first-order choice, we take §^ = [fy^] 00 , where ^ 



,(2) 100 



,(2) 



is the block-diagonal interaction giving rise to in (|2~77| ). Therefore 



(2) 



(3.32) 



where ^q 1 '* was determined in the previous section [formula ( |3.31| )1 and (fQ Z> is such that 



Vn° + VP 



V- (2) 



Inspection shows that many of the terms comprising V®° + Vj° are zero when the original 
perturbation Q is just a nearest-neighbor hopping term. In fact, <y? y is nonzero only for the 
following sets Y (we continue to use the letter X for (unordered) pairs of nearest-neighbor 
sites) : 

(i) Y = X, a nearest-neighbor pair: 



(2) 

Vox 



1 



Px- 



(3.33) 



(ii) Y = X U X' , where X and X' are different nearest-neighbor pairs sharing a site: There 
are two contributions, 

(3.34) 



(2) ~(2) , _(2) 

foxux' — V^OXUX' + V^OXUX' • 



The first one is 



~(2) 



J2 P< xux> [ad^ix 4 (ad5ix 3 (ad5ix 2 (Qx 1 ))) 

c 6 



P< 





XUX' 



(3.35) 
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where Cq is the set of sequences (X\, X 2 , X 3 , X 4 ) satisfying one of the following conditions 





= x 2 


= X; 


x 3 


= x 4 


= X' 


< x x 


= x 3 


= X; 


x 2 


= X A 


= X' 


. Xx 


= x 4 


= X; 


x 2 


= x 3 


= X' 



u 



same but 
X < — > X' 



(3.36) 



The second contribution is 



-(2) 

^OXVJX' 



o 4 



p 





XUX'i 



(3.37) 



&dS 2 (x i ux 3 ) (<32(Xiux 2 ) > 
with C 4 being the set of sequences (Xi, X 2 , X 3 , X 4 ) satisfying one of the following conditions 

U { 7 } ^ (3.38) 



X\ — X 3 — X ; X 2 — X4 — X' 
X\ = X 4 = X ; X 2 = X 3 = X' 



(Hi) Y — X U X' U X" U X'" forming a unit square of the lattice: For concreteness, let us 
assume that, in the sequence X,X',X",X"', each set shares a site with the next one, and 
X'" shares a site with X. Again, there are two contributions 



(2) _ ~(2) 1— (2) 

{ PqXVJX'\JX"VJX'" — { PoXUX>UX"UX>" ' V 9 0XUX'UX"UX"' • 



(3.39) 



~(2) 1 
^0IUX'UX"UX"' 



- V P° 

o r XUX'UX"UX"' 



C24 



The first term is the four-bond version of Q3.35| ): 

(adSWadSWQxJ)) Pxux'ux"ux"' 

(3.40) 

C 24 = {(A^A^A^X,) : X 1 ,X 2) X 3) X 4 e{X,X / ) X /, ) X w }}. (3.41) 
The second term is the four-bond version of (|3.37p : 



where 



1 



)0 



^0XUX'UX''UX"' — 2 ^ r XUX'UX"UX"' 

c 4 



adS 1 . 



2(X 4 UX 3 ) 



ux 2 ) 



pO 



(3.42) 



where C 4 is the set of sequences (Xi, X 2 , X 3 , X 4 ) satisfying one of the conditions 
[Xi = X, X 2 = X', X 3 = X", X 4 = X'") (J {3 cyclic permutations}. 



(3.43) 



Each of the operators defining the terms (|3.33|) -( |3.42|) can be expressed in terms of a 
product of creation and annihilation operators. The expressions thus obtained can be sim- 
plified by using the spin operators defined in fl3.20| ). To avoid unenlighteningly complicated 
formulas, we write the final result only in the direction-independent case, t\ = t b ± = t± . 
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^From ( |3.31|) it follows that if t±b = t ± ^ = t±, and for Y = (xy), 



^0{x,y} 



2{t 2 + +t 2 _) ( q3c3 P { X ,y} \ 4t+t_ c±c± 
JJ \°x°y A / * TT °i D «i 



4 1 ' [/ 

whereas from ( |3.33| )- (|3.42| ) we obtain the following expressions: 
') For Y = (xy) , 

<p?Ly = w + + 6 ^ 2 -) - %) - ^ + s * s v 

We omit local projectors, because of property ( |3.23| ), and we denote 

C-L Q-L cl Cl I C"2 c2 / c+ Q- j_ C— C+"\ 

°x °y ■— °x D y °x J y ~ 9 ^ x V x V ' ' 



(3.44) 



(3.45) 



(3.46) 



(%%') For Y = X U X', with X = (xy) and X' = (yz): 



<£0{x,y,z} 



and 



fo{x,y,z} 



pO 

{x,y,z} 



U 3 



+ 



v2(t*+tt) , ^t\tl 
4 (f.t- + t 3 _t 



U 3 



( q3 c3 1 c3 q3 q3 q3 ' 



f/ 3 



,y.2 ^2 



pO 



pO 

{z,y,z} 



2 {t% + ti 



f/ 3 



_ In 8tit 2 
■* 4^ ' U 3 



(3.47) 



/ q3 o3 ± \ 1 ut '+ t '— / q3 q3 1 q3 q3 q3 q3 
'^x^ J J "i ^3 \ D x D y ~r D y D z D x°z 



*t 2 +t 2 _ ± ± A(t 3 + t_+t 3 _ t +; t(1 cfl c± , 

£/3 a z ' £/3 



pO 



(3.48) 



(m'j For Y = X U X' U X" U X'", wift X = (xy), X' = (yz), X" = (zw) and X'" = (zw): 

-(2) nO f 8 (*+ + ^ 



^0 {x,y,2,Ul} 



po r 

{a;,3/,z,u;} \ 



u 3 



q3 q3 q3 q3 
O x O y O z O w 



2 (t + + t_ ) / 3 Q 3 „3 Q 3 ^3 „3 Q 3 Q 3 ^3 Q 3 ^3 q3 _ ^ \ 
^3" "r T 'Jz'Ja, T >~> W D X O x D z O y D w ~) 



— (t^jt- + t 3 _t + ) (SyS^ S^S^J + (S^S^) SyS z + (S^Sy ) ^ ^ 4" ( ^ ^ ^ ) 



/ Q-L C-L j_ C-L C-'- _i_ C-*- C-'- _i_ Q-^- Q-^- 
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\tit 2 



- - — (S+Sy SzS w + S;,. SyS z 

4 (Sx^z) Sy$l + 4 (SyS^) SlS z - S^Sj- - SyS^ j P^^} 



f/ 3 
Atltl 



U 3 



(3.49) 



and 



w (2) P° 

ro{x,y,z,w} 1 {x,y,z,w} 



{ 



r4 {t% + t 4 



c/ 3 



( q c3 c3 c3 c3 c3 o3 q3 Q 



32 i 
+ ^3 + t-t+) (S z S w ) S x S y + (S x S y ) S Z S W + (S^ 5 2 ) S X S W + (jSj. S^) S^iy 



32 1 2 t 2 

H — (S x S y S Z S W + 5*+^ S+) 



u 3 

8t 2 + t 2 _ 
U 3 



8 {Sx^z ) SySfv + 8 (SySw) Sl^z + + SyS^ | -P^^^^} • (3.50) 



Two particularly interesting limits are the following ones. 

The symmetric model: Setting t + = £_ = t in the above expressions for and grouping 
operators which have the same support, we arrive at the expression for the effective Hamil- 
tonian of the standard one-band Hubbard model in the low energy sector, to fourth order in 
the hopping amplitudes: 



+ (~Jj ~~ IT 3 ") 2( S * ' Sy 4^) + ip ^ P {x,y,z} (Sx-S z ~ -) P{ x , y ,z} 

(xy) ' (xyz) 

|/ ' - 

— ~Tj3 ^2 P{x,y,z,w} \Sx ' Sj, + • S 2 + s z ■ s w + s w ■ S x 

1 x y J 

+ Sx ■ S Z + Sy ■ S W — -^j P{ X ^y tZtW } 
[(Sx-Sy) (S^.S^) + (S^.S^) (Sy.S^ (Sj;.S 2 ) (Sj,.S t 

U fwz\ 

\ xy J 



(3.51) 
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where the symbol (xyz) denotes that x and y are nearest neighbor sites and so are y and 
z, while 

{Ty} 

stands for four sites around a plaquette. Expression ( |3.51| ) has been derived 
previously. It has been reported, for instance in 0, |39|] (where the hopping parameter is 



allowed to depend on the direction), in @ (together with higher-order corrections, up to 



eighth order) and in | 31| . Our contribution amounts to showing that if one systematically 
uses local projections, unlike in the first three references, the remainder can be controlled 
rigorously. This is the new result, which follows from our results in ||. 

The Falicov- Kimball limit: Setting t + = 0, t_ = t, one obtains the effective interaction for 
the Falicov-Kimball model to fourth order. 



(2)^00 

fkJ 



f2t 2 2t 4 N ^ / 3 3 P{ X>y} 



2t 



\^ P° ( Q 3 Q 3 i C"3 c3 9 c-3 r-3\ p 

jj 3 r {x,y,z} \ D x D y i* D y°z z D x°z J r \ 



{xyz) 





{x,y,z} 



2 1 

\ p0 I q3 q3 i q3 q3 i q3 q3 , q3 q3 , o3 q3 , q3 q3 \ p0 

tt 3 y x y ' y z ' z w ' w x ' a? z y wj {x,y. 

U fwz\ 
\ x y J 



z,w} 



40 t 4 ^ o o o o 

\ q3 q3 c3 c3 
~T7T- 2.^ ^x^y^z^w ■ 



U 3 



(3.52) 



( w z\ 
1 IJ J 



For the purposes of determining (classical) ground states we can ignore the projection oper- 
ators, in which case the Hamiltonian becomes 



- EN 



^p0 

2 x 



2? 

IT 



18t 4 



At 4 fi f 4 40 1 4 

_l_ \ " c3 c3 , u ^ \ " c3 c3 j_ ^ u _ \ " c3 c3 c3 c 



3/1=2 



£/ 3 



|x-»| = V2 



t/ 3 



3 

tu ■ 



I is J 



(3.53) 



This effective classical Hamiltonian and its ground states, were studied in . In JTS 



a more 



general expression was obtained involving hopping amplitudes that depend on the direction: 
fox] _ -f.[y x \ _ i^Wjgi^^ _ —0 yx Such complex hopping amplitudes describe the influence 
of a magnetic field. Our expansion methods provide systematic control of remainder terms. 
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3.6 Higher-order perturbation theory 

While straightforward, in principle, the computation of higher-order transformed interactions 
is a tedious and error-prone process. (In fact, already the second order computations lead 
to monstrous expressions if hopping is allowed to depend on the direction as well as spin.) 
The use of computer codes to perform symbolic algebra becomes mandatory. To encourage 
industrious readers, we sketch some basic elements of the algebra needed to treat the resulting 
expressions. 

The different terms arising are products of two types of operators: (i) creation and 
destruction operators, and (ii) projections. We can therefore classify the relevant terms in 
three categories. 



3.6.1 Products of creation and destruction operators 

The algebraic manipulations of these operators are well known: one uses the anti-commut- 
ation relations ( |2.9| )-( PTTU| ) to obtain expressions in terms of the occupation numbers n xa = 
elixir an d the spin operators defined in fl3.22| ). Furthermore, differences of occupation 
numbers can be written in terms of if>J via ( |3.2 



3.6.2 Products of projections 

There is only one equation to consider: 

P'Pj = 1-6,,. i,j = 0,1. (3.54) 

A similar equation holds for projections on any set X. Another pair of potentially useful 
relations are 

P x p0 x = P °x > P xP~ = P 1 x , (3-55) 

whenever X G X. 

3.6.3 Products of projections and creation- and destruction operators 

This is the non-trivial part of the algebra. Perhaps, its most important relations are the 
intertwining relations 

c P° = P 1 c c P 1 = P° c 

^xa 1 x x ^xa ^xa 1 x x ^xa 

(3.56) 

r t po _ pi „t r t p 1 — po r t 

xa x x xa xa x x xa ' 

for arbitrary a — ±. They allow a purely algebraic verification of which terms are zero. 
We remark that projections of the form P x , with X containing more than one site, are not 
directly suited for the use of ( |3.56| ), because these projections are not products of single-site 
projections. We then use that 



P \ = i-p° = i- n*? 

x&X 



(3.57) 
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As an example, let us sketch the algebraic proof that 

[adS laX (Q° a } x ,)] 00 = (3.58) 

for X = (xy) and X' = (yz), with x ^ z; a fact only parenthetically mentioned after ( |3.19| ). 
Indeed, the left-hand side of ( ft.5£ ) is a combination of products of the form 



P{x,y,z} P{x,y} C ya C xu P{x,y} ^{y,z} C zcr> C V<r' ^{y,z} ^{x,y,z} ' (3.59) 

plus various x < — > y and y < — > z permutations. We first notice that the middle projections 
P\ can be removed. To see this we use (|3.57|) and then move, for instance, the resulting 
factors Py to the far left and far right through the intertwining relations ( 3.56[ ). Such factors 



emerge as Py and cancel with the corresponding factors Py at both ends [relation fl3.54 )]. 
We are left with 

P{x,y,z} C ya c xo C za , C ya i P{ x ^y^ z } , (3.60) 

where we have also used (|3.55| ). To see that such a term is zero, use again the intertwining 



relations to show, for instance, that as there is an odd number of creation and annihilation 
operators involving the site z, the leftmost P° can be written as P], on the right. It then 
annihilates with the corresponding P 2 °. 

There are many other useful relations involving products of projections and creation and 
destruction operators. Besides the previous formulas ( |3.23|) , (|3.26| ), (|3.27|) and (|3.28|) / (|3.29|) , 
we mention 

n xa P° x = [l-n^P* (3.61) 
n xa Pl = ^pi (3.62) 

and 

<L><LP£ = (3.63) 

c^c^Pl = (3.64) 

for any a, o' = ±. By the intertwining relations (|3 .56 ) the identities (|3.61|) -( ^.64|) can also 
be written with the projections on the right or on the left and the right. 

4 Low-temperature phase diagrams 

4.1 Foreword on quantum Pirogov-Sinai theory 

The purposes of Sections ^| and ^ was to construct a unitary operator, U^(\), with the 
property that the conjugated Hamiltonian 

flf(A) = U£\\)H A (\)U£\\)* 

= H$(\)+Vi n \\) (4.1) 
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of some quantum lattice system is "block-diagonal" to order |A| n+1 . By this we mean that 
the matrix elements of (A) between ground states, or approximate ground states (low- 
energy states) oiH^(X) are of order n + 1. Our method to construct uj^\X) has been (a 
convergent form of) analytic perturbation theory. 

In attempting to gain some insight into the structure of ground states and low-temperature 
equilibrium states of the Hamiltonian H^\X), it is tempting to argue that one should ne- 
glect the perturbation (A) and study the ground states and low-temperature equilibrium 
states of i?oA (A) — hoping that they are close to those of (A). Unfortunately, this hope 
is ujustified, in general. 

It is a well known fact that arbitrarily small perturbations can have a drastic effect on 
low-temperature phase diagrams. This is typically the case when the truncated part has an 
infinite degeneracy. For example, when there is a continuous symmetry, states are usually 
grouped in bands with gapless intraband excitations. An arbitrarily small additional term 
can split the lower bands and hence change the phase diagram. Such models are hard to 
treat rigorously; in particular, no such treatment is available for the (symmetric) Hubbard 
model, despite our control on the perturbation series of Section [|. The case in which the 
truncated interaction has an infinite degeneracy not connected with a continuous symmetry 
is usually more tractable. Often, entropy effects (quantum and classical) are the determining 
aspect, and there exist techniques to deal with them, as we will illustrate below. 

The simplest situation arises when the truncated interaction has a zero-temperature 
phase diagram involving only finite degeneracies. In favorable cases, we expect that each of 
these ground states gives rise to a corresponding ground state, and a low-temperature state, 
of the full, untruncated interaction. The detection and description of these phases, however, 
requires some additional machinery. Indeed, even at zero temperature, where only energy 
plays a role, we cannot rely on plain diagonalization procedures. First, there is no hope of 
achieving exact block diagonalization of the whole family of Hamiltonians {-£^a} needed to 
pass to the thermodynamic limit, barring miracles (that is, relatively trivial and uninteresting 
situations). And, second, the estimation of the effect of the "undiagonalized" remainders is a 
delicate matter. The remainder term, V^(A), makes a negligible contribution to the energy 
density. But it may have a decisive quantum entropic effect: It may alter the overall wave 
function by allowing "virtual quantum transitions" between classical configurations. The 
situation is even more involved at nonzero temperatures, due to the appearance of thermal 
fluctuations. 

Quantum Pirogov-Sinai theory 

Its bare-bones version [§, 0] applies to systems with a leading part that is "classical" , in 
the sense of being diagonal in a tensor product basis (Section |2.2| ). It is based on contour 
expansions. The Duhamel expansion is used to write partition functions and expectations as 
sums over (d+l)-dimensional piecewise cylindrical surfaces — the (d+l)-st axis corresponding 
to the inverse temperature (imaginary time) axis. Such "contours" have cylindrical pieces, 
corresponding to Peierls-like contours for the classical part of the interaction, and a finite 
number of section changes caused by the quantum perturbation. The objective of the theory 
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is to prove that, in the thermodynamic limit, these contours are sparse and far apart. This 
yields a precise mathematical representation of states as "classical ground states modified 
by quantum and thermal fluctuations" . In fact, both types of fluctuations are treated on the 
same footing. So quantum Pirogov-Sinai theory provides a unified description of ground- 
and low-temperature states. 

In brief, the theory establishes sufficient conditions for the energy to dominate over 
further quantum and thermal entropic effects. As a consequence, the latter can only cause 
small (and smooth) deformations of the zero-temperature phase diagram of the classical 
part. That is, the phase diagram remains stable under quantum perturbations and small 
raises of temperature. This does not mean that all phases remain stable throughout the 
phase diagram. For instance, small entropic effects can tilt the balance in favor of particular 
phase (s) of a coexistence region and cause the disappearance of the remaining ones when the 
thermal or quantum perturbation is switched on. But the theory assures that the unfavored 
phases reappear for slightly modified parameter values, that is, the perturbed system will 
also exhibit a coexistence manifold, albeit a little shifted and deformed. Moreover, the theory 
provides a criterion to determine which phases are stable, and, in fact, it presents a picture of 
what "stable phase" means. [As a matter of fact, below we shall resort only to this stability 
criterion. The part on the stability of phase diagrams will not be useful, as coexistence lines 
will, in all cases, involve infinite degeneracies.] 

As we discuss in more detail below, the conditions on the classical part required by the 
theory are of two types: First, it must lead to finite ground-state degeneracies within the 
region of interest, and, second, it must satisfy the Peierls condition which roughly requires 
the energy cost of an excitation to be proportional to the area of its boundary. The first 
requirement forces us to stay away from regions of infinite degeneracy. So, in principle, we 
would be unable to deal with the line /i + = /x_ of the phase diagram of Figure [I]. We shall 
see below that the quantum hopping actually helps us there. Similarly, models with low-cost 
energy excitations, like the balanced model, are out of reach. In fact, in [[12] it is shown how 
quantum hopping brings this model within reach. 

When applicable, the theory does a complete job: It accounts for entropic effects and 
provides a full description of states. Nevertheless, the theory has a limited scope In fact, it 
seems to be applicable only when nothing interesting happens as a result of the perturbation. 
The hidden card is our perturbation technique. In many instances, it produces a partially 
diagonalized interaction with a leading classical part incorporating crucial quantum effects. 
If such a leading part satisfies the right hypotheses, the subsequent application of Pirogov- 
Sinai theory will prove the existence of phases exhibiting truly quantum-mechanical features. 
In fact, a slight extension of Pirogov-Sinai is needed in these cases, showing that transitions 
within the excited band have a negligible effect, regardless of the order. Such an extension, 
based on a more careful definition of contours that distinguishes "high-to-high" transitions 
from the rest, is presented in ||. Below, we prove that hopping produces quantum sym- 
metry breaking among the ground states of the on-site Hubbard interaction, which brings 
the degeneracy down to a finite one. This is an instance of "quantum entropic selection". 



Another example, presented in ||12|| , is the "quantum restoration of the Peierls condition" 



taking place in the balanced model. 
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Incidentally, the thermal entropic selection is a better known phenomenon, which, in 
some cases, can trigger a complicated pattern of phases with cascades or staircases of phase 
transitions. This effect has been the object of extensions of Pirogov-Sinai theory based on 
the notion of restricted ensembles ]3], |3^, [19], f|, . The extension of this theory to quantum 



systems is a promising direction of research. The work in can be considered as a step in 
this direction, though it does not treat genuinely degenerate quantum interactions and only 
resorts to restricted ensembles as an alternative way to exhibit degeneracy breaking. 

Quantum Pirogov-Sinai theory acts, thus, as the link needed to pass from diagonalization 
—at the level of interactions — to phases. It is a powerful tool that can be used as a black 
box: Once the system is seen to be thermodynamic within its range of applicability there is 
nothing else to do; the theory provides a flexible and complete description of the resulting 
phases. We know of no other comparably general and versatile approach to the study of 
quantum statistical mechanical phase diagrams. We see the methodological difference be- 
tween quantum mechanics and quantum statistical mechanics. In the former, the ultimate 
goal is to diagonalize, exactly, the Hamiltonian of some system. In the latter, an exact diag- 
onalization is impossible. The goal is then to do the minimal perturbative diagonalization 
needed to bring the system within the reach of some theory enabling us to determine the 
phase diagram, like our Pirogov-Sinai approach. Once such a theory takes over, we are done: 
All the required information is at our disposal; any further diagonalization is a waste of time. 

4.2 Summary and hypotheses 
4.2.1 Scope of the theory 

Quantum Pirogov-Sinai theory applies to interactions of the form (|2.18 ), that is, 



$ = $ cl + Q , (4.2) 

where $ cl = i s a finite range classical interaction (classical in the sense of Section |2.2[) , 

and Q = {Qx} is a possibly quantum interaction. The operators comprising both parts 
depend on a finite family of perturbation parameters A = {A 1; . . . , A^}. For our applications 
below, $ cl will be the leading part <&o^ of a transformed interaction of an appropriate order 
n. The theory has two consequences: 

(PS1) Stability of states: Under suitable hypotheses, it determines the different periodic states 
of the full interaction, for small values of the parameters A and small temperatures. 
Each of these states can be traced, in fact, to some ground state of <3> cl , from which it 
differs in the addition of a diluted gas of thermal and quantum excitations. In such a 
case, the corresponding ground state of $ cl is said to be stable at the given temperature 
and values of A. 

(PS2) Stability of phase diagrams: It establishes (sufficient) conditions under which the zero- 
temperature phase diagram of the classical part remains "stable" under the addition of 
the quantum part, Q, and/or increasing temperature. That is, the manifolds where dif- 
ferent phases coexist are, for small temperatures and quantum perturbations, smooth 
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and small deformations of the corresponding coexistence manifolds of ground states of 
the classical part. 

The theory requires two sets of hypotheses: (i) the classical leading part must have a 
finite ground-state degeneracy and satisfy the Peierls condition explained below, and (ii) the 
matrix elements of the quantum part must be sufficiently small. 

4.2.2 Hypotheses on the classical part 

Let us discuss the hypotheses on the classical part first. The finite degeneracy refers to 
the presence of only a finite number of periodic ground states for $ cl . In the present set- 
ting (m-potentials with minimal bond energies normalized to zero) these are configurations 
. . . ,u p such that § c y(uji) = for all bonds X. The Peierls condition refers to the well 



known generalized Peierls contours |$5], ^4], |37|]. They are constructed by means of sampling 
plaquettes 

Wj a) := {y E Tiy : \ Vi - x t \ < a, for 1 < i < v} , for x E . (4.3) 

The set of sampling plaquettes where a configuration u does not coincide with any of the 
ground states Ui is called the defect set of u>. The contours of u are pairs 7 = (M, u M ), where 
M — the support of 7 — is a maximally connected (with respect to intersections) component 
of the defect set. The radius a of the sampling plaquettes must be larger than the range of 
$ cl and the period of each of the u>i, so that, knowing the set of contours, one can univocally 
reconstruct the configuration oj. In particular, when $ cl is a transformed interactions $o™' ) , a 
must be larger than the radius R of the plaquettes used to define protection zones (Section 



23). 

For models whose interactions and excitations are determined by nearest-neighbor con- 



ditions, like the Ising model, one can use the "thin" Peierls contours |33], [10], [L3[] traced 



with segments midway between pairs of adjacent sites. But in more general cases one must 
resort to definitions as above to ensure a one-to-one correspondence between configurations 
and families of contours. The energy of a configuration is, in the present setting, the sum 
of energies, £(7), of contours 7. Here, £(7) is the energy of the configuration having 7 
as its only contour. A model satisfies the Peierls condition if each contour has an energy 
proportional to the cardinality of its support, 3(7). That is, 

£(7) > ks(7) (4.4) 

for some k > called the Peierls constant. 

In this paper, we consider a more detailed Peierls condition based on the distinction of 
two levels of excitations — low- and high- lying — of the classical part. The distinction is 
made through a suitably defined family of local projections Py- With respect to this family, 
a configuration u exhibits a high-lying excitation in Y if P y uj 7^ uo. For the applications of 
this paper, these projections are exactly those defined in Section [2.2| , because the low-lying 
excitations of $ cl = are among the local ground states of the original $o- The high- 
energy defect set of a configuration uj is the union of the plaquettes where it exhibits 
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a high-lying excitation. A system satisfies a two-level Peierls condition if for each contour 7 
we have 

E(i) > Ks(7) + Ds(7 Mgh ) (4.5) 

where 7 hlgh is the part of the high-energy defect set contained in 7. This condition means, 
in particular, that high-energy excitations lead to an additive excess energy measured by 
D > 0. 

Given that $ cl is allowed to depend on the perturbation parameters A, such a dependence 
is also expected for the Peierls constants k and D. 

4.2.3 Hypotheses on the quantum part 

First, we decompose the operators of the quantum interaction, Q = {Q x }, into "low^low" 
(££), "low— >high" (£h), "high— > low" (h£) and "high— >high" (hh) components: 



Qx '■- p x Qx P% 

:= (1 X -P° X )Q X P X (46) 

Qx := P°xQx(lx-P x ) 

Qf := (1 X -P X )Q X (1 X -P° X ). 

The conditions on the quantum part are that there exist sufficiently small numbers e ai 
and S such that, for a, 7 = £, h 

HQxl < e a ,6^ x K (4.7) 
The region of validity of the two-level quantum Pirogov theory depends on the parameter 



Quantum Pirogov-Sinai theory converges provided n is sufficiently small. 

The original (one-level) quantum Pirogov-Sinai theory discussed in 0, Efl corresponds to 
taking D = e M = e ih = e hh = 0. 

Interactions obtained as a result of the partial block-diagonalization procedure of Section 
0, contain terms proportional to powers of X/D; (A measures the strength of the quantum 
perturbation, and D is proportional to a typical energy-denominator). Condition ( |4.7|) thus 
imposses the bound 

max(A ,4) < 8 , (4.9) 



D' 

where 5 is as in ( f4.7|) ; (5 must be chosen small enough for the perturbation expansions to 
converge). The size of the h£-matrix elements of the quantum perturbation, {Q x = Q x (A)}, 
in the Hamiltonian H^\x) is estimated by 

\n+l 

WQxW - WQxW - -^r^ x) - n ' 1 - (4-10) 
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Comparing (|4.10|) to (|4.7|) we find that 

\n+l 

^ = = 0(^ TT ) • (4.11] 

As the "high— >high" part of the original interaction survives unaltered, 

•As 



O(-) . (4.12) 



4.3 Stability of phases 
4.3.1 The basic criterion 

Our criterion for the stability of phases is based on the use of cluster-expansion technology 
to construct some objects /{,..., f' k — where f[ is associated to the ground state uii of $ cl — 
called truncated free energy densities. They are defined provided that 

e := maxfe^ ^ , rj) < £ , (4.13) 



for some (small) constant £o that depends on parameters like the range of <3> cl , dimensionality 
of the lattice, size of the sampling plaquette, and dimension of the on-site Hilbert space. 
Within the convergence region, the truncated free energies are analytic functions of j3 and 
of any parameter on which the interaction has an analytic dependence. They determine the 
stability of phases in the following way. 



Stability criterion Assume that A is such that $ (A) has a finitely degenerate ground 
state and satisfies the Peierls condition (j.5), and assume that Q(X) satisfies hypotheses 
( JTp . Then if, for some values of f3 and A within the region ( |4.13| ), 



Re/p(/x) = mm Re /J (/i) , 

for some p with 1 < p < k, the following holds: 

(i) /p(/i) coincides with the true free-energy density of the system; 

(ii) the infinite-volume limit 

tr? Ae"K P 



lim 

A/-z d tr^ e-^ H A 
exists for any local operator A; 



(iii) 



(A)7 A - (uj p \A\u p ) 



< HAIIlAlOfe) 



for any operator A G Ax, {X is a finite subset of the lattice). 



(4-14) 



(4.15) 



(4.16) 
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The notation H A P stands for the Hamiltonian with "external condition" u> p , that is, the sum 
of all $x with X fl A ^ 0, but allowing only matrix elements between vectors that coincide 
with Up outside A. The symbol tr^ p indicates a trace over the space of such vectors. 

Note that, in (|4.15| ), H^ p is the "effective" Hamiltonian, rather than the original Hamil- 



tonian defining the model. The effective Hamiltonian and the original one are unitarily 
equivalent, the unitary conjugation being given by an operator ujy\X) as in fl2.81|) , ( |2.84|) , 



( |2.85| ). In order to select a stable phase, we impose boundary conditions u p outside boxes A 
on the effective Hamiltonians rather than the original Hamiltonians. One could reconstruct 
"boundary conditions" for the original Hamiltonians corresponding to the boundary condi- 
tions that we impose on the effective Hamiltonians. (Typically, "boundary conditions" for 
the original Hamiltonian will constrain configurations not only outside A but also insided A 
but near dA. We shall not attempt to provide details concerning the map from boundary 
conditions for the original Hamiltonians to boundary conditions for the effective Hamiltoni- 
ans. Boundary conditions do not have an operational meaning here (corresponding, e.g. to 
certain experimental conditions), but are mathematical devices to select stable phases that 
should be chosen in as convenient a way as possible. In fact, this is the conventional role of 
boundary conditions in statistical physics. When we describe a magnetic material in terms 
of a quantum Heisenberg model we are working with an effective Hamiltonian; (there are no 
explicit exchange interactions present in the fundamental Schrodinger Hamiltonian!). But 
we do not hesitate to impose boundary conditions directly on the Heisenberg Hamiltonian, 
rather than on the fundamental Schrodinger Hamiltonian of the material, in order to select 
stable pure phases (e.g. the direction of spontaneous magnetization). 
Enough! Let's get down to business! 

We shall say that there is a stable Wp-phase whenever ( [4.14|) is satisfied. 

We shall not enter here into the technicalities of the actual proof of this criterion, which 
can be found in ||. Instead, we present some general comments on the basic construction. 

The above criterion is proven by resorting to a low-temperature expansion obtained by 
iteration of Duhamel's formula for the exponential of the sum of non-commuting matrices. 
This adds an extra continuous variable ranging from to (3. There is one such expansion 
for each ground state of $ , involving a "sum" of terms labelled by the sites of the lattice 
("spatial variables") and the continuous "inverse temperature" or "time" variable. Every 
term can be labelled by a piecewise cylindrical surface contained in 2 d x [0, 0\, which can be 
decomposed into connected components called quantum contours. Outside the surface, the 
space is filled with the corresponding ground state of $ cl , and the contours are transition 
regions whose interiors are occupied by different ground states of $ cl . The contours have 
spatial sections formed by the usual classical contours described above, which grow cylindri- 
cally in the "time" direction until there is a sudden change of section due to the action of a 
quantum bond. The hypotheses on the quantum part ensure that the weight of this contour 
decays exponentially with its area: The exponential cost of each "vertical" cylindrical piece 
is given by the Peierls condition Q4.5p , while each "horizontal" change of section is penal- 
ized by a term ( |4.7| ) according to the type of transition involved. Each contour has a finite 
number of section changes. 
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We must distinguish between "short" and "long" contours. The long contours traverse 
the whole "time" axis from zero to (3. Therefore, their cost is at least exp(— const (3), and 
they disappear in the limit j3 — > oo. In particular, a purely classical system would have only 
contours of this type which, in fact, would be straight cylinders with no section changes. 
"Short" contours appear and/or disappear at intermediate values of the continuous integra- 
tion variable. Thus they involve some changes of sections, and their cost is proportional to 
5. Such contours survive the limit (3 — > oo and contribute to the quantum ground state of 
the full interaction $. 

We see that, in this approach, thermal and quantum effects are put on a similar foot- 
ing. They are both sources of entropy associated to different contour geometries. A low- 
temperature equilibrium state can be visualized ctS cl SGcl ' configured of the corresponding 
ground state of $ cl plus the fluctuations represented by contours. In particular, the short 
contours can be interpreted as "ground state fluctuations" . At low temperatures and small 
values of 5, these fluctuations are dilute, because their large energy cost overwhelms the 
"entropy gain". Thus expectations in such a state differ little from expectations in the 
associated ground state of $ cl . 

4.3.2 Stability and symmetries 

Two types of symmetries can be distinguished. First, there are symmetries that leave each 
term of the interaction (or of an equivalent form of it) invariant. Such a symmetry is 
associated, for each finite volume A, to a unitary operator Sa such that 

SaH^SX 1 = H s K Uv . (4.17) 

[The vector Suo p is easy to visualize at the level of classical configurations. In a quantum 
statistical mechanical formalism, it is defined, for instance, via limits of St}uj p \ Examples of 
these symmetries are the spin-flip symmetry of Ising models in zero field, or the particle-hole 
symmetry of the Hubbard model at half filling. 

On the other hand, there are symmetries associated to operations on the lattice, like 
translations and rotations. They map each term of the interaction into a different term, 
but leave the whole interaction invariant. Each such symmetry is defined by a bijection 
T : Z d — > Z d that yields natural bijective maps between configurations — (Tu) x = 
and between local operators — (T$)x = ^t-^x (for notational simplicity we denote the 
different transformations with the same symbol). The latter map does not leave Hamiltonians 
invariant, but connects Hamiltonians corresponding to translated volumes: 

TH jp T -i = Ht-I k ■ (4.18) 

In general, a symmetry, R, can be a composition, ST, of symmetries of the previous two 
types. 

The criterion of stability of phases presented above respects symmetries. Indeed, we say 
that a symmetry R connects a family of configurations {u>i, . . . , Uf} if 



R l ui = for 1 < i < £ - 1 



(4.19) 
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Then, under the hypotheses of the preceding subsection, we have: 

Ground states of $ cl connected by a symmetry of $ are 
either all stable or all unstable. 



(4.20) 



We briefly sketch a proof of this fact, as it is not explicitly given in || or in ||. The key 
object to look at is the partition function H P (V) for piecewise cylindrical, bounded regions, 
V, of 7L d x [0,(3]. This object is defined, through the contour expansion, as the sum of all 
allowed contour configurations inside V compatible with the boundary condition uj p . For 
regions V of the form A x [0, 0\, with A a finite region of Z d , one has 



2p(A) 



tr A p e 



(4.21) 



Looking carefully into the definition of H(V), it is not hard to conclude that its behavior 
with respect to symmetries is similar to that of ( 4.21|) . Namely, if S a symmetry of the first 
type, 

Sa*(V) = Z„ p (V) , (4.22) 
and, if T is a symmetry of the second type 



-TV 



(4.23) 



Here T _1 V refers to the space-time region V transformed by the map T _1 . 

These partition functions are involved in the following characterization of stability: A 
ground state u p of $ cl is stable in the region ( |4.13| ) if, and only if, there exists some constant 
C such that 

~-' ,n ' < exp(C\dV\) , (4.24) 



for every bounded region V C TL d x [0,(3] and 1 < j < k. Here \dV\ indicates the area of 
the external boundary of V. 

The criterion ( [4.20|) is now easy to verify: Assume that some ground states U\, . . . ,aj£ 
are connected by a symmetry. If the symmetry is of the first type, then by ( |4.22| ) 



E e (V) 



(4.25) 



for all piecewise cylindrical V C %. d x [0, (3\. Hence (|4.24j) is either verified by all p — 1, . . . , £ 
or by none. 

Consider then a symmetry of the second type, and assume, for concreteness, that u>e is 



stable. Then (|4.24j) is verified for p = £ and all j, in particular for j 
(l<q<£) [see ($T%\. But then, from QOg ) we have 



q where uj q = T ^ 9 - ) 



E^T-i^V) 



< exp(C|d^|) . 



(4.26) 



The last inequality follows from the assumed stability of uj^ and the fact that |<9(T\^)| = \dV\. 
Therefore, for any 1 < q < £, and any 1 < j ' < k, 



z»t(Y) 



< exp(2C|^|) 



(4.27) 
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where the inequality follows from the stability of uji [eq. ( [4.24 )1 an d from ( |4.26|) . We conclude 



that all the ground states U\, . . . ,u>£ are stable (with an exponential constant 2C). 
4.4 Stability of phase diagrams 

More generally, one studies families of interactions $ M parametrized by a finite set of pa- 
rameters ji = {fix, . . . ,/i s } — typically fields or chemical potentials — and one is interested 
in determining the corresponding phase diagram, that is, in obtaining a catalogue of the 
different phases present for different values of /i. If for the values of /i under consideration 
the system satisfies the hypotheses described above, one can apply quantum Pirogov-Sinai 
theory and determine the phase diagram on the basis of condition ( |4.14 ), now involving 



/^-dependent truncated free energies f . The cluster expansion tells us that the difference 
between these truncated free energies and the corresponding energy densities is of the order 
of the parameter Ep\ introduced in ( |4.13| ). Therefore, given the smoothness properties of /', 



one would expect that the phase diagram, for (3 and A in the region (|4.13|) , differs little from 
the phase diagram of $ cl at zero temperature. 

This can, indeed, be proven, under some minor additional hypotheses: Chiefly, (i) bounds 
similar to ( |4.7] ) but involving the partial derivatives dQ^x/dfii, and (ii) a hypothesis of 
regularity of the phase diagram of The latter roughly means that the parameters fi 
completely break degeneracies among the ground states, so the Gibbs phase rule is satisfied. 
The detailed hypotheses can be found in || Section 5.2]. The conclusion is: 

If the hypotheses sketched above are satisfied for an open region O of the space of parameters 
fi, then the phase diagram for f3 and X in a region of the form is regular and is a 



smooth deformation of the zero-temperature phase diagram of in O. The displacement 
of the different coexistence manifolds is of the order of ep\. 



5 Phase diagram for the Falicov- Kimball regime 

We now apply the phase-diagram technology described in the previous section to asymmetric 
Hubbard models, whose perturbation expansion has been discussed in Section |3|. Let us 
denote 



sup 



sup 



(5.1) 
(5.2) 



(with a slight abuse of notation). Let us see how further features of the phase diagram 
appear as we consider higher orders in the perturbation expansion. 
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5.1 Phase diagram to order t° 

To this order, the classical leading interaction correspond to a (formal) Hamiltonian 

H ci(o) = uJ2n x+ n x _-Li + J2n x+ -Li-J2 n *-- ( 5 - 3 ) 



We see from Figure [1] that it trivially satisfies the hypotheses of Section |4.2.2| — in fact for 



the one-level theory — as long as one avoids the lines where the different regions with a 
unique ground state intersect. In the latter we loose the required finite degeneracy. The 
conclusion is, therefore, that different ground states remain stable for chemical potentials in 
the open regions of uniqueness. The range ([4.13 ) shrinks to zero as we approach a coexistence 



manifold, because so do the different Peierls constant (which are of the order of the energy 
of the lowest excitation). This phase stability remains true upon the addition of electron 
hopping or any quantum perturbation which is translation-invariant and satisfies rj < Eq (for 
instance small ion- hopping) . 

5.2 Phase diagram to order t 2 /U 

To this order we can exhibit more features of the phase diagram within the shaded region in 
Figure |2|. The transformed interaction was obtained in Section |3.4| [Formula ([3.31 )1. It nas 
a leading classical part 



o 



# cl « = U^n x+ n^ + ^\t [ -^{S^-^)-h^Sl, (5.4) 

x (xy) x 

whose periodic ground-state configurations are given in Figure |j. This interaction satisfies a 
two-level Peierls condition. The low-lying excitations have one particle per site and involve 
a Peierls constant 

k ~ At 2 /U ±h + 0(t 4 /U 3 ) . (5.5) 

The excess energy for configurations in the non-singly occupied band gives an extra Peierls 
rate 

D = m&x(U — no , Ho) — 0(i7) (5.6) 

(Ho defines the boundaries of the shaded region in Figure Q for concreteness we are consid- 
ering Ho ~ U). We see that the leading part (|5.4|) satisfies the hypotheses required for the 



classical part in Section [4.2.2| as long as 

At 2 /U±h + 0{t 4 /U 3 ) > 0, (5.7) 

that is, except within bands of width 0(t A /U 3 ) centered at the lines h = ±4t 2 /U. Moreover, 
the summability of the transformed interaction is guaranteed only if D ^> k. This condition 
determines the limits of the shaded region in Figure |2|. 

The "low^low" transitions must include at least one ionic jump, hence a factor t + . 



Therefore [see eg. the quantum correction in ( |3.31 )] 
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where 5 has been chosen sufficiently small to guarantee the convergence of the partially 
diagonalized interaction. From this, ( |4.11| ) (with n = 1) and ( |4.12| ) we see that the quantum 



perturbation satisfies hypotheses ( }4.13|) if 



y,^<?o, (5.9) 

where Eq is a small number that shrinks to zero when k tends to zero (and at the boundaries 
of the shaded region in Figure ||) . 

We conclude that quantum Pirogov-Sinai theory [together with the symmetry criterion 
( |4.20| )1 proves that for nonzero but small values of t, the phase diagram of the model is as 
in Figure The ground states correspond to the configurations shown in the figure plus 
quantum fluctuations (=short space-time contours). Furthermore, they remain stable for 
small temperatures (long contours appear sparingly) and small values of t + (new contours 
appear, rarely, involving the new type of quantum transitions) . This picture is valid outside 
an excluded 0(t 4 /[/ 3 )-vicinity of the lines h = ±2t 2 /U where the Peierls constant k cannot 
be guaranteed to be positive, and inside the shaded region of Figure ^ where This 
stability persists under the addition of arbitrary quantum perturbations satisfying ( |4.7| ) that 
do not break translation invar iance. 



5.3 Phase diagram up to order t 4 /U 3 

In order to analyze in more detail the excluded regions of widths 0(t 4 /?7 3 ), we need to 
consider the next order of perturbation theory. For the direction- independent case (t b _ = t), 
the leading classical interaction, to this order, is given by the terms in ( |3.53|) . At this point, 
the hardest part of the analysis is, by far, the determination of the ground states of this 
complicated classical part. Fortunately, this has been done before, see [|5|, [Tj|]: the ground 
states correspond to the configurations depicted in Figure ||. We see that the classical part, to 
fourth order, satisfies the hypotheses required by Pirogov-Sinai theory, except in the vicinity 
of four lines of infinite degeneracy. In regions at a distance 0(t 6 /U 5 ) from these lines, the 
Peierls constant k is of order t 4 /U 3 . The gap D is still given by ( |5.6| ). 

Moreover, the quantum part has coefficients Em and given by (^4.11|) with n = 2, and 
Ehh given by ( (4.12[ ), while 

^°(tw) < 51 °> 

[see eg. the quantum corrections in ( |3.45| )- (|3.50|) 1. 



As a consequence, the stability criterion of Section [4.3. 1| is applicable, for large (3 and t 
satisfying ( |5.9| ) with a smaller e . Together with the symmetry considerations ( |4.20| ), this 
implies the stability of the ground states of Figure [| — except at the excluded bands of 
width 0(t 6 /U 5 ) around the coexistence lines — under the addition of temperature, ionic 
hopping [within the limits impossed by ( |5.9|) 1, and any other translation- and rotation- 
invariant quantum perturbation satisfying ( f4.13p . The phase diagram of Figure |5| is thus 
obtained. 
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Figure 8: Set B x formed by the bond (xy) and its protection zone for the 3-band Hubbard 
model 



A similar analysis can be performed for cases in which the hopping depends on the 
direction. In particular, in the presence of magnetic flux, ie. if t_ = texp(i6 yx ) with 
Qyx £ (0, 2tt), our block-diagonalization procedure yields, to order t 4 /U 3 , an effective classical 
interaction with terms given in [18, Formula (3.11)]. As shown in this reference, the phase 
diagram for the interaction with flux involves the same ground states as in the flux-less case, 
but with deformed manifolds of infinite degeneracy. The Pirogov-Sinai approach proves 
the stability of these ground states, under quantum and thermal perturbations, except in 
excluded regions of width t 6 /U 5 around the lines of infinite degeneracy. 



6 Perturbation expansion for the 3— band Hubbard 
model 

In all the models considered in the previous section, the Coulomb interactions between the 
particles were strictly on-site. In this section we apply our perturbation method to the 
3-band Hubbard model, as an example of how to proceed when the leading interaction has 
nonzero range. 



6.1 The original interaction 



The Hamiltonian Hqa in (|1.2| ) is of unit range, since it consists of both on-site and nearest 
neighbor interactions. This necessitates the use of "protection zones" [see Section [2.2|| for 
the proper treatment of the local character of the operators that arise in our perturbation 
expansion. We define i?-plaquettes, W x , through fl2.23|) with R = 1, and sets B x through 



( (2.24Q . In particular, for X := (xy) , x G A,y G B, \x — y\ = 1, B x is given by the set of 
sites shown in Figure |8] (or a rotation and/or reflection of it). 

It consists of two copper sites and four oxygen sites. We restrict our attention to param- 
eter values in the following ranges: 



U d > e d ; A > ; U d > U p > U pd > ; U d > A ; t pd > ; U d > t pd , 



(6.1) 
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and normalize the ground state energy of Hq\ to zero. The unperturbed Hamiltonian is 
given in terms of an interaction <3> — {^oy}) nonzero, if Y — (xy) denotes a pair of nearest 
neighbor sites such that x <E A and y G B. The interaction $ x is defined as follows. 



OY 



U d n° 



Up TT'y-^ fly 

U pd n d x nP 



n d x _ + £d (n d x - 1) liY = {x} , x e A 
HY = {y},yeB 
if Y = (xy) , x e A,y e B 



+ e p n p y 



(6.2) 



For parameter values satisfying ( |6.1D , the energy of a configuration on Y, with respect to 
the interaction $oy, is minimum (and equal to zero) when there is a single hole at x and 
none at y. Hence in a ground state configuration of Poa, each oxygen site is empty while 
each copper site is singly occupied by a hole, the total number of holes in A being equal to 
the total number of copper sites, \A\. However, the ground state of H \ has a 2' A '-fold spin 
degeneracy. 

The interaction corresponding to the 3-band Hubbard Hamiltonian H\ ( |1.2j ), with the 
normalization introduced above, is given by $ = where 



$y = $ y + Q Y , 

and Qy = unless Y is a pair of nearest-neighbor sites, in which case 

Q{xy} = t pd [p\ a d XCJ + d\ a p X(T \. 



(6.3) 



(6.4) 



We define suitable projection operators on the Hilbert space Ha as in Section |2.2| . Let Pg 
denote an operator which projects onto (local) ground states of the interaction $ - The 
projection operators Pg and P| are defined by ( |2.25|) and (|2.26| ) of Section f2~2"[ . Using the 
resulting partition of unity 

1 = P° Bx +P 1 Bx +Pb x i (6-5) 

we can decompose the perturbation interaction Q x of fl6.4| ) as in ( |2.37| ). It is clear from the 
structure of the lattice [Figure [J and the definition (|6.4|) of Qx that = 0. Hence 



Qx — Q°b x + Qb 



(6.6) 



where Q° B l x and — defined in ( p.39| ) and ( |2.40| ) — are linear in the hopping amplitude 
t pd . Also, due to our choice of ground-state energy normalization, 



oo 

0B> 



(6.7) 



6.2 First-order perturbation for the 3-band Hubbard model 

As explained in Section |2|, we first search for a unitary transformation U^(t p d) = exrj(t p dSi] 
with 



1B, 



(6.8) 
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which eliminates the first-order off-diagonal terms Q°g of the perturbation interaction Q 



From ( 2.5U| ) we have that 



where 
with 
Hence 



E = and E 1 = U pd + A, (6.10) 
A:=e p -e d . (6.11) 

Sib x = jj d + ^ { P b x Qb x p b x - p b x Qb x p b x } (6-12) 

for X = (xy), x G A and y G B. A straightforward calculation along the lines sketched 
in Section |2.3| shows that the transformed classical interaction has nonzero terms only for 
nearest-neighbor pairs X. These terms are [recall (|6.7| )1: 



$ ol = -P Bx ^S lBx {Q^ x )P° Bx 



1 

2" 

+2 

(C/ pd + A) X L ^ a xaPy<rPya a x(7 J I! x 

^CE(»i(i-»y)^ 



In the last line we used the identity 

P°b x ( E "1,(1 - ng„))Pg x = ( E = L ( 6 - 14 ) 

0-= + ,— (7= + ,— 

We conclude that, to this order, the change in the classical interaction amounts to an 
irrelevant shift in the (local) ground-state energy that does not introduce any new effect. In 
particular, it fails to reduce the spin degeneracy. We need to go to the next order of our 
perturbation scheme to find non-trivial contributions. 

6.3 Second-order perturbation for the 3-band Hubbard model 



^From Section \LM we obtain that, in second-order perturbation theory, the leading classical 
part, $J, 2) = [^o\ 00 , is of the form 

4 2) = + tp® , (6.15) 
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where the nonzero terms of the latter are 

1 



(2) 
V0B X 



p0 



adS 1Bx (adSiB x (adS 1Bx (Q Bx ))) 



(6.16) 



for X a nearest-neighbor pair, and 



(2) 

VOBy 



J2P°b y [ad5 1Bx4 (ad5 1Bx3 (ad5 1Bx2 (g i xi ))) 



1 



P 



By 







adS 1 ' 



2(S X4 USx 3 ) 



UB X2 ) 



y > 



(6.17) 



for K = XUX', where X = (xy), X' = (yz) are pairs of nearest neighbor sites with x,z £ A 
and y £ B. The assignments C6 and C4 in the last two sums are those defined in ( |3.36|) and 
( p.38| ) , respectively. 



The one-bond contribution (|6.16|) is again an uninteresting energy shift: 



(2) 



(U pd + A) 



- P° 

3 B{x,y} ' 



(6.18) 



The two-bond contribution ( 6.17[ ) can be cast in a more familiar-looking form by making use 
of spin operators at copper sites, [see (|3.20|) 1: 



(6.19) 



One obtains 



(2) 
VOB 



(U p d+ A)< 



Jr. 



with 



J, 



cff 



{x,y,z} 



4& 



(s«.s, 



{%,y,z} ' 



1 

Pd 



2A + £L 



(6.20) 



(6.21) 



(P prf + A) 2 

^From ( 6.15| ), ( |6.13| ), ( |6.18 ) and (|6.20| ) we recover the known fact that the effective Hamil- 
tonian for the 3-band Hubbard model in the low-energy sector, to order 4 in the hopping 
amplitude t p d, is given by a S = 1/2 antiferromagnetic Heisenberg model on the square lat- 
tice of copper sites. This is the same interaction obtained for the one-band Hubbard model 
to order t 2 [third line in (|3.51 )]. This observation justifies, in part, the use of the simpler 
one-band model in studies of magnetic properties of undoped cuprates. Our contribution is, 
once again, to be able to provide convergent estimates on remainder terms in the transformed 
interaction. 
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